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SYNOPSIS 


TH30EETIGAL STUDY OP THE SDSOTRONIC STRUCTURE OP GOVAIENT 
SOLIDS AND DOPED SEMICONDUCTIONS 

by 

PRADSSP KUMAR KHO¥ASH 
Department of Physics 
Indian Institute of Technology 
Kanpur 

The aim of the present work is to study theoretically 
the electronic structure of a number of ionic, covalent and 
semiconductor solids, using MS Ko;SCP method. The valence 
electron energy levels and the equilibrium charge distribution 
of all these systems are analyzed to get informations on 
their bonding characteristics, optical properties and 
electrical conductivity. 

The MS Xj^method is briefly described in Chapter I. 

The method considers a cluster of atoms or ions s copped out 
of the solid and replaces the effect of the rest of the solid 
by putting proper charges on the boundary of the cluster 
(i.e., on the Watson sphere). This is expected to simulate 
the localized physics in the solid to a good approximation. 

The cluster is divided into atomic regions, interatomic 
region and extramole cular region. One assumes a muffin-tin 
approximation to the potential for each atom in the cluster. 
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The potential is spherically averaged within each atomic 
region and in the region outside the Watson sphere. In 
the interatomic region the potential is volume averaged. 

The exchange potential is taken in the local X^. approximation 
(p'^/5) Slater. The basis wave functions as solutions of 
Schrodinger equation are obtained in the atomic regions, the 
interatomic region and in the extramole cular region. The 
Wronskian on the boundaries are set to zero in order to 
obtain the energy eigenvalues, and consequently the eigen- 
functions. The calculation is perfornBd self-consistently. 

The name multiple scattering theory comes from the fact that 
the wave function in the interatomic region is built with 
the waves scattered from the different atomic sites. The 
exchange parameter a are due to Schwartz. 

In Chapter II, the MS method of calculation is 

applied to systems like ionic solids (NaCl, liF, HaP), 
perovskites (M = Pe, Ni)), rutiles (MP 2 (M=Pe, Co,Ki,Zn) ) , 

metallic oxides (1%0, I&O, VO) and alloys (Pe— Te). Along with the 
energies the charges within each atomic region are calculated. 

A charge transfer diagram is defined to understand the chemical 
bonding (ionic and/or covalent) in these systems. The 3d 
degenerate level of the cation breaks up into t^g and 
levels in the ligand field. The crystal field splitting, 

10 Dq, is determined for crystals with 3d cation and compared 
with experiment wherever data are available. The transitions 
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between two of the valence levels are calculated wherever the 
corresponding XPS or optical spectroscopic experimental data 

are available. Finally, the Mb'ssbauer isomer shifts are 
calculated for iron containing samples. 

The self-consistent charge distribution study in ionic 
solids like NaCl, NaP, LiP shows a transfer of charge from 
the central ion to the ligands along with some transfer to 
the interatomic region. This indicate that these well known 
ionic compounds are not perfectly ionic but show a marginal 
covalent character. Por NaCl the calculated charge distri- 
butions are in very good agreement with the experimental 
results. The energy levels of liP and NaP bunch into two 
groups with predominantly P 2s and P 2p character. The 
energy difference of the P 2s - like and P 2p - lilce states 
in liiP is 1.43 Sy as compared to the experimental result of 
1.54 Ry obtained by XPS. 

The energy levels and the charge distributions are, 
then, calculated for fluorides and oxides of the different 
transition metals. The agreement of the energy difference of 
1.36 Ry of the P 2s - like and P 2p - like states in PeP 2 with 
the XPS result 1 .54 Ry is encouraging. The crystal field 
splitting, 10 Dq, is also calculated for the transition metal 
fluorides. The Mg 2s and % 2p levels a,re identified and 
their energy difference 2.85 Ry is found to be in fair 
agreement with the XPS data 2.38 Ry. 
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The Mo'ssbauer isomei' shifts in all the iron containing 
systems, namely KPeP^, PePg g alloy are calculated 

and compared with the existing results from the experiments. 

The first step tov;ards the analysis of the covalent 
nature of the systems is to calculate the covalency parameter 
from the coefficient of the ligand term in the extended 
eigenfunction of the uppermost valence level. These values 
are matched with the obtained from the analysis of tte 
neutron diffraction and nuclear magnetic resonance experimental 
data in terms of the molecular orbital theory. Satisfactory 
agreements are obtained for EPeP^ and KNiP^. Por KPeP^ 
our results are = 5.58% as against the neutron diffraction 
data = 6 . 24 %. Por KNiP^ they are 7 . 55 % and 5 . 11 % 
respectively. The value of for llgO is 6.67 percent which 
is in good agreement with the SKDOR value of 6,50 percent. 

The concept of pseudoatoms is then brought in and 
the charge transfer from and/or to these atoms are calculated 
to plot the charge transfer diagrams for all the solids 
studied by us. With the help of these dia.grams, the measures 
of ionic bonding strength and covalent bonding strength for 
these systems are defined. One finds that the natural state 
of the non-net all ic solid is a mixed ionic- covalent state. 

In NaCl the ionic is the dominant bonding and covalent bonding 
is marginal whereas in PePg the covalent bond strength is 
very large and ionic bonding is marginal. The transition metal 



perovskites and rutiles generally fall within these two 
limits. However, they are predominantly covalent. The 
oxides of the transition metals also show very large values 
of covalent bond strength with different degrees of minor 
ionic character. is an exception to this general character 

in that it has a small charge transfer from the anion to the 
cation. This type of charge transfer, referred to as inverse 
ionic character, is typical of semiconductors. MgO is known to 
behave as a semiconductor. Inverse ionic character along 
with large covalency is seen in all semiconductors e.g.Si, 

InP, GaAs and transition metal impurities in InP. 

In Chapter III, the MS X^SCP method was applied to 
semiconductors, especially transition metal doped InP, with 
a view to study the localization of the levels within the 
band gap. The cluster model has been applied to the canonical 
semiconductors such as Si and GaAs. It is well known that 
Gr make GaAs semi-insulating. A similar role played by 
Pe in InP is expected and so the self consistent field energy 
levels and their eigenfunctions are worked out. Detailed 
calculations predict a highly localized impurity level for 
Pe in InP. 

Por the other transition metal impurities in InP, 
the impurity levels are either closer to the va.lence band 
or conduction band edges. These states thus hybridize with 
the host and become delocalized as compared to Pe in InP. 

Very recent experiments on Pe and Go in InP confirms a higher 



localization of the former leading to a higher resistivity. 

The calculations of the different charge states of 

DBH 

Pe in InP are also performed. The tg -level, arising due 
to broken bonds, find its place in the mid gap for neutral 
Pe , whereas it moves closer to the band edges for the other 
two charged states. The charge transfer is minimum for 
neutral Pe. The isomer shift for different charge states are 
calculated but unfortunately no experimental data are 
available to compare with. However, these results also confirm 
that neutral Pe makes InP semi- insulating. 

The electronic structure of the alkali metals as 
impurities in InP are also calculated. It is seen that 
they interact weakly with the host thus pinning the impurity 
level in the gap. 

It is predicted that GaAs may soon supercede Si in 
technological importance. Since the devices made out of 
InP are much faster than those of Si, it is expected that 
InP is going to be a keen competitor to GaAs as a device 
material. This work forms the theoretical basis for the 
analysis of the conductivity properties of transition metal 
doped InP. 

In conclusion, we hope that the present work has 
extended our knowledge of the electronic structure of covalent 
solids and semiconductors with impurities, especially of their 
bonding characteristics. 



,..for an atom of F electrons, the solution is a function 
of 3N variables, and even if it were possible to evaluate 
such a solution, to any degree of numerical accuracy required, 
no satisfactory way of presenting the results, either in 
tabular or graphical form, is known. It has been said that 
the tabulation of one variable requires a page, of tv/o 
varia.bles a volune , and of three va.riables a library; but 
the full specification of a single wave function of neutral 
le is a function of seventy-eight variables. It would be 
rather crude to restrict to ten the number of values of each 

variable at which to tabulate this function, but even so, 

78 

full tabul8.tion of it would require 10 entries, and even 
if this number could be reduced somewhat from considerations 
of symmetry, there would still not be enough atoms in the 
whole solar system to provide the material for printing 
such a table. 

D.R. Hartree, Reports on Progress in Physics, Vol.11 

(London; The Ph^rsical Society, 1947), p. 113. 
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MSXflSCF METHOD 

IIKRODUCTION 


In the present -work we have used the Multip;j.e 
Scattering Self Consistent Field (MSX SCP) method [1 ] 
to study theoretically the electronic structure of covalent 
solids and semiconductors with transition metal impurities. 

The valence electron energy levels and their eigenfunctions 
have been obtained and analyzed to get informations on their 
bonding characteristics, optical properties and electrical 
conduct ivity . 

The MS Xq;SGP method considers a cluster scooped out of 
the solid. The effect of rest of the solid is simulated by 
enclosing the cluster by a sphere, called the Watson sphere, 
and putting proper charges on the boundary of the sphere to 
make the cluster electrically neutral. The potential within 
each non- overlap ping atomic sphere is taken in the muffin-tin 
form. The potential is spherically averaged within each 
atomic region and in the region outside the Watson sphere 
and volume averaged in tte inter atomic region. The exchange 

1 /3 

potential is taken as P as given by Slater [2], The basis 
wave functions for a solution of the Schrodinger equation 
are obtained in the atomic regions, interatomic region and 
in the extramole cular region. The Wronskian at the boundaries 
of the region are then set to zero in order to obtain the 
energy eigen values. The calculation is performed self- 
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consistantly. 

We applied this method of calculation to systems like 
perovskites (KMP^ (M = Pe, -Ni)), rutiles (MPg (M=Pe , Go, Ni, Zn) , 
ionic solids (NaCl,liP, NaP), metallic oxides ( I%0 , MnO , VO ) 
and alloys (Pe-Te). Along with the SOP energies the charges 
within each atomic region are also calculated. A charge 
transfer diagram was defined to get an understanding of the 
chemical bonding (ionic and/or covalent) in these samples. 

The 3d degenerate level breaks up into t^ and e levels 
in the presence of a crystal field. The splittings are found 
to be in good agreement with the experiments. The isomer 
shift was calculated for tte iron containing samples. 

The transition metals are known to introduce levels 
deep in the band gap of the host semiconductors [3-5 ]» called 
the deep impurity states. They play a very important role in 
modifying the electronic properties of the semiconductor. 

The III-V materials have a very promising technological future 
since the devices made out of them are much faster than those 
made out of silicon. We perforned a MS /^SCP calculation for 
pure and doped InP for the first time. The transition metals 
used as impurities were Or, Mn, Pe, Co, Ni and Cu. The 

electronic structure shows that the dangling bond hybrid 

DBH 

level, t 2 , is pulled down as the atomic number increases, 

DBH 

owing to the greater attractive potential. This tg -level 
is exactly in the middle of the band gap for Pe in InP. This 
provides a clue to the semi- insulating behaviour of Pe in InP, 
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A very recent experiment (1984) shows the resistivity of .Pe 
doped InP to he larger than Co doped InP. The different charge 
states of Pe in InP was studied to have a detailed understanding 
of the semi- insulating behaviour. The alkali-metals were 
also doped in InP to find out if they passivate the dangling 
bonds arising due to vacancy. 

1,1 One Electron Approximation tethods 

In a system of N electrons and k nuclei (assumed to be 
at rest), each electron moves in the electric field produced 
by all (N-1 ) electrons and all the nuclei. For K greater than 
a very small number, a rigorous solution is impossible. By 
averaging the density of (N-1 ) electrons over their complex 
motions and finding the electric field from this average 
density helps one in reducing the many body system to a single 
body problem. The individual electrons then move in a fixed 
external field and their wave functions as found from the 
solutions of SchrSdinger equation (SS) are called orbitals - 
atomic orbitals for atomic case and molecular orbitals for 
molecules or crystals. 

The atomic orbital problem is somevrhat simple because 
of its spherical symmetry. However, the SS for even one 
of the MO'S is difficult to solve. The only basic simpli- 
fying feature there is the point group symmetry of the 
molecule or the cluster which simulates the basic features 
of the crystal. Consider a system of nuclei with coordinates 
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and charge and masses and electrons xvith coordinates 

r.. The time independent SE in atomic units (e.u., ’^=e=m =1 ) 
1 6 

is 


H(R^,Z^,M^,rini.(E^,9l) =BupB^,r. 


( 1 . 1 ) 


We restrict to spin free non-relativistic Hamiltonian. 
Our hasic interest is in the electronic orbitals. Using Born- 
Oppenheimer approximation 


H .U , 
el el 








'A' 


( 1 . 2 ) 


where 





S 2 
1 A 



+ 


2 2 J_ +22 

i>-ir. . A>BR 


(1.3) 


The electronic wave functions depends 

directly on the electronic coordinates and parametrically 
on the nuclear coordinates. Throughout we will assume that 
there are no ’slow' electrons or ’fast’ nuclei which might 
lead to strong electron-phonon coupling. 


The quantum chemical methods for the solution of 
equation (1.2) can be broadly divided into four classes. 


1 . Semi Empirical ^thods; 

, Model Hamiltonian 
. Adjustable molecular parameters 



6 


e.g., . Hflckel - Simple tight binding 

Extended HU eke 1 

. CNDO 

2, Ab-initio methods; 

. Choose a form of U (one or many determinants) 
, Choose a set of basis function 
. No further fundamental approximation 
e.g., . LCAO-MO-SCI’ (Hartree-Hock) 

. Configuration interaction (01) 

, Generalized valence bond (G?B) 


3. Psudopotent ial nethods; 

. Like ab-initio (or local density) but core electrons 
are replaced by pseudopotentials 
e.g., . SCP, Cl, GVL, 

4. Local density methods; 

, Approximate local exchange correlation potential 
e.g., . Xa, LSD 

, Roots in density functional theory and theory of 
electron gas 

, Various ways of solving the electron SS 
e.g., . LCAO, DVLOAO, SW method. 

Instead of considering the whole crystal we can focus 
our attention on a small cluster whose repetition in space 
2 ?eproduGes the crystal as a whole. The symmetry properties 
of the cluster help in simplifying the complexity of the 
problem. The most effective part of the cluster calculation 
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is that it can s<juee^e out the localized physics quite 
accurately by approximating the effect of the rest of the 
solid in a suitable manner. 

The importance of cluster calculation relating 
different fields [6-H] ai'e shoxvn in figure 1.1. 

The various Xq; LSD are tabulated as 

Xp, -SW - Scattered wave 

Xq; “MS - iviultiple scattering 

- Muffin-tin potential 

- Numerical solutions in spheres 

- Partial waves 

- Most rapid method 

- Accurate orbitals and spectroscopic 
quant it ie s 

“ Does not yield total energy curves 

Xq. -DV - Discrete variational method 

ICAO (Slater orbitals) 

1 /3 

- fit P and P ’vjith auxilliary functions 

- Numerical sampling for matrix elecents 

2 

- Quite rap id -work goes up as m (ia=Fo. of basis sets) 

- Very large number of sampling points 
needed for total energy 

ICAO Xq - ICAO (Gaussian orbitals) 

1 Z’? 

- fit P and P ' xvith auxilliary-fjmsrtions 


Analytical integrals 
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Fig. 1.1 

Some fielc^ of oppUcation 6f cluster science. 
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3 

- Slowest of Xg, methods - work goes up as m 

- Accurate total energies 


1.2 Hart re e and Hart ree -Pock Methods 

n 

The many electron wave funct-i^n was written down by 
Hartree as product of the individual wave functions 

N 

i U > =. n 5 (1 .4) 

i=1 ^ 

The one electron Hartree equation is given by [15] 


[- + V^(1) + \(1)] U^(1) = Vi(1) 


wnere 


V,(1) = 


2Z 


+ 


U (2) U (2) 
S n. / ^ 


3 0 


1 


1:2 I 


dx, 


and, 


9,(1) = 


U.(2) U,(2) 

"i ^ Ti; ""2 


(1 .5) 

(1 . 6 ) 

(1 . 7 ) 


A 

V is the field arising from the nuclear charge and 
the total charge density of the electrons in the system and 
V, Gorifects the coulombic term for the self interaction of 
the electron. 


The Hartree model ignores the spatial correlation in the 
motion of any pair of electrons produced by their instan- 
taneous coulomb Impulsion. In Hartree-Pock (HP) the coulomb' 
correlations are again ignored but the spin correlation is 
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introduced through determinantal wave functions in place of 
equation (1.4). This is called statistical correlation as 
the electrons obey Fermi -Dirac statistics. This correlation to 
the coulomb potential made to get an acceptable SGF is called 
exchange correlation potential which contains those spin 
dependent exchange effects that are characteristic of anti- 
symmetry of the wave function and lead to magnetic properties 
of atomic systems. The neglect of coulomb correlation is 
corrected in a large measure by the inclusion of statistical 
correlation. 

The Hartree-Pock one electron equation [16,17] is 
- r + V ("I ) 1 U. (1 ) - 2 n.7U*(2) U.(2) , ^ 

Uj_(l) (1.8) 

The last term on the left hand side now includes in 
addition to self correlation term 



(1 .9) 


(1 . 10 ) 


characterized as exchange term. The only non vanishing 

contribution in equation (1.10) arises from the summation over 

the orbitals U. which has the same spin as U . , owing to the 
3 ^ 
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orthogonality constraints. 


V. 


This term can be written as 
- ? n. / U. (1 ) U. (2) 


HPX^ 


(1) = 


ir-,- Tr 


U.(1 ) U. (2) di:. 
0 1 ^ 


* / 

U. (1 


U^(1) 


( 1 . 11 ) 


This exchange correlation potential is the potential 

of a charge whose magnitude is one electronic charge removed 

from a hole surrounding the location of electron 1 . This 

cnarge which we shall call the "exchange correlation charge", 

t h 

has the same spin as that of the i — spin orbital and has a 

charge density at the position of the first electron which is 

just great enough to neutralize the total electronic charge 

t h 

of the spin of the i spin-orbital at that point [18]. This 
hole termed as Term! hole [ 2 } accounts for the correlation 
of the like spin electrons only. So the last term on the left 
hand side of equation (1 .8) should be corrected for the 
Coulombic correlations of unlike spin elnctrons. 


1 Z Pree Electron Exchange Approximation 

The Hartree-Pock exchange potential for different 
occupied orbitals are similar in many respects. So the 
idea is to form an universal exchange potential by suitably 
averaging these individue.l exchange potentials. 

The physical nature of average exchange potential avoids 
calculation of exchange integrals altogether. 



Previously the electronic structures were calculated 
by Thomas-Permi model where the charge density was treated 
statistically. Dirac refined the method by adding the 
exchange terms proportional to the cube root of the charge 
density. The resulting Thomas -Permi-Dirac (TPD) model 
gave a simple good approximation to the potential in an 
atom. However, it failed to describe binding between atoms 
in molecules and solids. In practice TPD model was never 
used as more accurate HP equations could be solved. The 
major defect of tfce TPD model was the crude treatment of 
the kinetic energy (PE) of the electron which is proportional 
to the 2/3 power of the charge density. Slater in 1951 came 
up with another model where the ES was treated as in HP model 
but the exchange term was approximated statistically. This 
corrects the drawbacks of TPD model and is simpler than HP 
model. 


In a non uniform system. Slater proposed the form of 
exchange potential to be 


ysxch(^) ^ -3[^ p(r)] 


1/3 


(1 . 12 ) 


The inhomogeneous electron gas problem has been solved 
using this local electron correlation potential [2,19-21], 

A rigorous theoretical justification for such local density 
potential methods is provided by the theorem of Hobenberg 
and Hohn [22 ]. According to them the ground state energy of 
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an inhomogeneous electron system is a unique functional 
of density P(r). Por a given external potential V{r) 

E[P(r)]= / V(r)P(r) d^r + E[P(r)] (1.13) 

This theorem merely proves the existence of a unique 
functional and offers no insight into it’s exact nature [ 23 ]. 

A number of functionals have been proposed to suit various 
types of many electron systems [24,25 ]. Equation (1.13) can 
be written in detail as 


E[ P(r)] =:P V(r) P (r) d^r + ^ d^r’ 

+ Tg[ P(r)] + P (r) ] (1,14) 


VThere the second, third and fourth term on the right 
correspond to the classical coulomb interaction, EE of 
the non-interacting many electron system with density P (r) 
and the electron correlation termi, A formal similarity 
between the Hartree equation and equation (1 .14) leads to 
a set of one particle equations |26] . 


[_ ^ 
L 2m 


V(r) + V^Jr) ] U. (r) = e. U. (r) 

(1.15) 


With, 


\Pr) = 


3E r (r) ] 

XC ^ / J 


3[P(r)] 


(1 .16) 
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All the many body effects ar& included in V which 

xc ’ 

unlike Hartree-Fock potential is local. 

As a first approximation the many electron system may 
be regarded as a homogeneous electron gas. The exchange 
correlation potential is then the 'exchange only' potential 
of the homogeneous electron gas [2o,21 ,27 ] . The exchange 
energy of the electron gas 


S 

xc 


[P(r) ]= - ^ f ( 


J/3 


8 ft 


[ P(r)] 


4/3 


d^r 


(1 .17) 


gives the exchange potential 

]= -2 [tV (1.18) 

comparing equation (1,12) (Slater’s value), and equation (1.18) 
(Gasper, Kohn and Sham (GKS) value [26,28]), one finds the 
difference of a factor 3/2. The difference is due to the 
fact that in GKS method the electron gas exchange is taken 
first and then the variation is made while in Slater's 
method [2 ] the variation is taken before the approximation js 
made. So finally a parameter a was introduced such that the 
final form of was obtained as 

= -3« [(3/47^ P (r)]'’/5 (1.19) 

a= 1 for Slater and 3/2 for Gasper, Kohn and Sham approxi- 


mat ion. 
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1.4 Xq; Method 

The approximation of Hartree-Fock method with Slater's 
exchange (HFS) is called the Xq, method. (X = Exchange and cc 
is a parameter introduced in the exchange term). method 

has gained a wide range of applicability in predicting 
localized properties through cluster calculation. The earlier 
approximations were restricted to atoms only. Pratt [29,30], 
Wood and Pratt [31] x\forked on transition meta.l atoms. Results 
of some non-self consistent calculations were found to be in 
good agreement with experiments but some of the self consistent 
results differred greatly. The method was then applied to 
molecular systems. Since 1965 much of the energy band 
calculations have been to determine tte values of a , which 
is unity, 2/3 or something in between [32]. The numerical 
method for a system which had neither the rotational 
symmetry of an atom nor the periodic symmetry of a crystalline 
solid was developed by Johnson [33-56 ] , The Bloch theorem 
is missing but some feature of the E3CB model of energy 
band calculation exist. 
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13 The Choice of a 


The X(y_ total energy depends on the exchange parameter 
Unfortunately, the total energy E is not a variational 
with respect to a . In fact, to a good approximat ion,E 

^ec 

is a linear function of c with large negative slope. If 
are the self consistent solutions of the X^, equations 


then, 


E [U^“ ] - <T(U^')rH! T(u2') > 


where H is the exact Hamiltonian and T (U^^ ) is the wave 
function derived from . This expression will be a 
variational with respect to 


3 S [ U^ ] 


:«= a 


Schwartz [37 ] tabulated the values of oc for different 
atoms and found that its value varied from §,7 for light 
atoms to 0.69 for heavy atoms except for hydrogen ( a = 0.9) 

A more practical way of optimising a is to calculate 
the KE and PE separately for the atom and choose a such 
that the Virial ratio is equal to unity [38 ] 


•2V [Uj_ ] 

__ 

I [u. ] 
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The defined in this manner shows a similar 

variation to ^ and the two values are almost equal for 

most of the atoms. This means that it is possible to choose 
an a such that Egp(oc) is close to HE limit and is 

close to 1 . 


By choosing a close to the HE limit and 

close to unity, we find the PE and EE close to the HE limit. 

Xa 

Since these terms are quite sensitive to the changes in , 
this optimization makes X(v and HE orbitals almost 
identical except for the virtual orbitals. 


A third way of optimizing ais to set E (o:'^)= E^^. 

% JlJ hJ:) 

In general ct should be diffeient for different orbitals. 
Schwartz has shown that the values of a vary only slightly 
as one goes from one atomic configuration to another (as in 
shifting an electron from the s to the d shell in transition 
element, or in removing an outer electron). This leads one to 
hope that it will not be too inaccurate to use the sa.me 
value of a for an isolated atom and for the same atom 
incorporated in the crystal. By setting different a for 
different orbitals no significa.nt effect in the model was 
observed. Things become complicated and the Virial theorem 


is not satisfied. 
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1.6 Method of Calculation 

The method is based on the division of complex or solid 
into polyatomic clusters. Each cluster contains the repre- 
sentei.tive basis atoms of the complex or the solid. 

The cluster containing same or different types of 
atoms is enclosed in a bigger sphere called the Matson sphere. 
The atomic sites form region I, the space in between the 
spheres form region II and outside the Watson sphere is 
termed as region III (Figure 1.2). 

To start witn, the intra-ionic or atomic potential 
within each ion or atom of the cluster with 110 mesh points, 
iS' interpolated to find the potential at 44l point mesh. 

Both the ionic potential and the eigen values corresponding 
to each of the atomic orbitals used as input for interpolation 
were taken from Hermann and Skillman [ 39 ]. 

The cluster potential at any point r is expressed 
as 

¥(r) = ? (Ir - R.I ) (1 .20) 

3 J 

where ¥• are the free atom SCF-]^potent ials centered at 
d 

positions R.. The potential within each atomic region I 
can be expanded about its own center in terms of spherical 
harmonics 


Vj (r) = I \{r) I^(r) 


(1 .21 ) 
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F!g. 1.2' . , 

Division of o molculof cluster into {D otcmic, 

{11} Inter otornic , and CMI) ext rcmolecubr regions* 



.20 


where L = The radial component can be 

evaluated by applying a , theorem by Lowdin [44 ]. Wfe then 
apply the muffin-tin approximation. According to this 
approximation we retain only the spherica.lly symmetric 
term 1=0 (1 = 0, m = O) in expansion (1.21 ) in each atomic 
region and treat higher order terms as perturbation. Hence 
the model potential is 


Y^(r) = V (r. ) 

o^ 



= 0 


b. 

1 


(1 . 22 ) 


where b^ is the sphere radius, includes contributions 

of all the other atoms to the site i. Thus the effect of 
the overlapping charge densities is partially included in 
spherical model. Electron-electron coulomb correlation is 
not within our model unless we deal with the total wave- 
f unction. The potential in the interatomic region is 
volume averaged 




V(r) dr 


(1.23) 


Vfe are concerned with the solution of Schr’ddinger 
e quation 


[ 



+ V(r)] T 


(r) 


Sf (r) 


(1.24) 


written in Hydberg units in region I, II and III. Within 
each spherical atomic region I of radius b-, we expand the 
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wave function in single centre form 

( 2 :) = 1 0<r<b^ (1.25) 

wiiere c£ are the partial wave coefficients to be determined, 
and Yjj( 21) real spherical harmonics. The functions 
(EJr) are solutions of radial equation 

[ _ d ^2 d_ ^ iXlill + v^(r) - E] R^(E|r) = 0 

r r 

( 1 . 26 ) 

for spherically averaged potential 'V^(r). The solutions are 
generated by outward numerical integration. 

In the region outside the fetson sphere, the solution 
is represented as 

'^'^^^(r) = I c°R°(E;r) ^^^(r) °° (1.27) 

where b^ is the radius of outer sphere. E° (E;r) are radial 
wavef unctions for outer spherically averaged potential, and 
are generated by inward numerical integration. 

In the intersphere region, the potential energy is 
volume averaged and Schr'ddinger equation is 

(V^ - E - Yjj) (r) = 0 (1.28) 

For E < its solution in the multicentre partial wave 
representation is given by 



22 


(r ) - 

ilV£; 




O . 


(1 .29) 


where i^(x) and modified spherical Bessel function 

and modified spherical Hankel function, respectively 


X= (1.30) 

and 

Si = I “ (1.31) 

J J 

Similarly for E > solutions of equation (1.26) can be 

written in terms of ordinary spherical Neum-ann and Bessel 
functions. We shall consider only E < case. 

In equation (1.29), the first term may be thought of 
as outgoing spherical wave scattered by regions of spherically 
averaged potential energy. The second term may be interpreted 
as incoming spherical wave scattered by extramole cular 
region. 


ivfe can also write the intersphere wave function as 

^ (1*32) 


where, 


II 


(r) = I (Xr.) Y^(rp 


(1.33) 


is the incident wave on atom j and 

®II (r> = £ 4 4 (Xrp Yj,(rp 


(1.34) 
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is the wave scattered from the j — atom. 

Thus the wave incident on atom j can be related to 
the waves scattered by all other atoms and interatomic 
potentials 2^ ^ ^ to the v/ave scattered by extramole cular 
potential. Thus we can express the coefficients in 
terms of the coefficients and Bor this, we male 

use of expansion theorems like 


(XI £2 " 




2 I (1-L’ ) 


k^„ (Xr^ ) (r^ ) i^, (Xrg) Y^, (rg) 


ri rg 


(1.35) 


and 


4 (Xlrg -r^O Y^(r2 - r^ ) = 4^ (-1)^^^' Ip,(L|L') 

, I 

(X r ^ ) Y^„ ( r ^ ) k^ , (X Pg ) Y^ , ( rg ) 


< ^2 

( 1 . 36 ) 


wheire integrals 


I^„(L|I’) 


S'rt n 

f d 0 j' sine d© Y^„(e,^) Y,(e,i^) y^, (e,i 2 f) 

00 J • 


(1.37) 


are nonzero only for following conditions 


FI - 1’ I ^ 1" < 1 + 1' 


(1 .38) 
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1" + 1 + 1’ = even integer 


(1 .39) 


Substitution of expansions into (1.29) and comparison 
with equations (1 .32) - (1 .34) yields set of equations 


"L “ j' L' "Lit 


. (B) 3^, (E) 


( 1 . 40 ) 


where 


G^', (E) = (1- 6. .,)47t(-1)l+l’ 2, I^„(l|l') k^„(XE^^,) 




(1.41) 


SJ^,(E) = 4^(-1)-^'"'- l„ Ip,(L?l') i^„ (><R^o)Yp. (R.q) 


(1 .42) 
(1 .43) 

(1 .44) 


R. 

— 

R., - R. 



“O' -0 

R-^ 


R - R. 

-OO 


-0 


For the outer spherical boundary, the 'incident waves' and 
'Scattered wave' intersphere wave function can be written as 

/ (r) = I Bi 4 (Xi’o) \ (£o) (1-45) 


f (r) 
^ II 


L ^1 ^I^-o^* 


(1 .46 ) 
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Equation (I. 45 ) giyes outgoing waves with respect to origin, 
singular at origin (r^ = 0 ) and equation (1.46) represents 
incoming waves with respect to the origin. A similar compari- 
son as for different spheres gives 


B 


0 

1 


S 

j' 


s 

L’ 


oOd’ 


(s) 4! 


(1.47) 


Now we equate the logarithmic derivatives of the intersphere 
wave functions with the corresponding values inside the spheres 
at various sphere radii b. and b , which gives following 

J U 

relations 

4 = *1 (S) (1.48) 

4 = '"l (E) 4 (1.49) 

where , 

tf(E) 

and 

t°(E) 

The square brackets are used to denote Wronskian 

- B(x) 

dx ^ dx 


[k( (Xb ), R°(Ei b )] 


[i,(Xb^), e" (E; bp] 


(1.51 ) 


[1 (Xbp, rJ (E;bp ] 
[4 (Xb.), (E; b ) ] 


(1 .50) 


[A(x), B(x)] = A(x) 


( 1 . 52 ) 
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Using 

identities 

[ij_(Xb), (Xb) ] = 

Xb 

(1 .53) 

and 

[jl(Xb),Bj_ (Xb) ] = 

1 

— 2 

Xb 

(1 .54) 


with equation (1.48) and (1.49)y following additional relations 
are obtained between partial wave coe.ff icients ; 

= (-1)^*^^ Xb^ [ Ij^(Xbj), Bj(E;b.) ]C£ (1.55) 

= (-1)^+^ Xb^[E° (E;b„), 4 (Xb^)]C° (1.56) 


The scattering on a given site of the molecule must 
be compatible with the scattering at different sites. The 
set of compatibility relations for the partial wave coefficients 
are obtained by combining equation (1,48) and (1.50) with 
(1.40) and (1.47) respectively. We get a set of linear 
homogeneous equations 


2 s 




II’ 



2 

L’ 


^LL' 


(E) Aj, = 0 


(1.57) 


S 2 
O' 1' 


^li' H 


2 

I’ 


6 

IL' 


[ti (S)] 



= 0 


where , 


(1 .58) 
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rT"^ (S) 

LL< 







i I 


LL' 


(S) 


(1 .59) 


The secular matrix of the problem is thus 




(s) 


-sj°, (E) 




The matrix elements are real and symmetric because real 
spherical harmonics are used and conditions (1.38) and (1.39) 
are satisfied. 


If the cluster has a point group symmetry, instead 
of spherically symmetrized wavef unctions we construct basis 
sets of the irreducible representation of this point group. 
The secular matrix then breaks up into blocks and becomes 
simpler to solve. 


1 .7 Advantage and Disadvantages of Xq; Bfethod 

The main advantage is the simplicity of the model 
and computational speed. It simulates the self-consistent 
environment of an electron by a single local potential. 

The Xq. orbitals are as accurate as double zeta basis set 
of the first and second row atoms and probably better for 
larger atoms which involve electrons with 1> 2 . 
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The orbital independent potential defined by Xq, 
method leads to a better one electron excitations of a 
system. Both the occupied and unoccupied (virtual) orbitals 
of N electron system are under the influence of the potential 
which simulates the effect of other (F-1 ) electrons. The 

ground state virtual eigen values are usually good descri" 
ption of the one electron excitations. A better first order 
picture of the spectrum is obtained if a transition state 
calculation [41] is done where one-half the electron is 
removed from a valence orbital and promoted to a.n excited 
orbital. 

The method also satisfies Virial theorem, Fermi 
statistics and the Hellmann-Feynmann theorem independent of 
a . This is convenient for calculating the force on a 
nucleus directly in terms of a three diirensional integral 
rather than through the evaluation of the six dimensional 
integral necessary for the total energy. 

The Xq. model can be considered to be an orbital 
generator which provides a convenient basis to be used 
as a starting point on a more accurate many body approaches. 
The method is rapidly convergent. There are no multicentered 
integrals to be evaluated. It takes care of the orbital 
relaxations. The energy expansion leads to Permi-Dirac 
statistics. The total energy can be found as a function of a 
sterio chemical geometry, Por periodic solids, the cluster can 



29 


be choosen as an unit cell. Application of Bloch's theorem 
can lead to band theory. 

The disadvantage is that the total energy does not 
correspond rigorously to a total wave function. The connection 
is lost with the approximation in the exchange part. The 
second disadvantage is the presence of a , the exchange 
paranEter, whose choice is not obvious. There is also a 
numerical disadvantage in this model. The exchange potential 
is a nonlinear function of the charge density. There are 
two usual common misconceptions; 

(i) Self interaction is disregarded - which is not true, 
(ii) Treatment of a free electron gas of constant 

density - can be had from a simple dimensional 
argument, quite independent of the perfect gas 
model [42] . 

The application to linear molecules turned out to 
be poor and the overlapping sphere method was introduced [43]. 

Nevertheless, since the cluster model is not dependent 
on the assumption of long range order, problems such as 
bonding of impurities and defects in crystals and the localized 
electronic structure of covalent materials are within the 
scope of the theory. 
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In general, the method does obtain a good estimate 
of the wave function and makes the quantitative treatment 
of any molecule or nonmetallic solid tractable. 

1 .8 Basic Structure of Program and Inputs 
The program consists of four parts; 

(i) Generation of atomic potentials 

(ii) Calculation of starting molecular potentials 

(iii) A non-SGB program to calculate the starting 
energy levels 

(iv) Calculation of self consistent set of energy 
levels and normalised wave functions. 

tfe shall briefly describe the parameters and various 
subroutines used for each part, 

(i) Atomic potentials; 

A Hermann-Skillman program is used to compute 
the atomic orbitals for the atomiS with a slight 
modification to incorporate for different atoms 
as given by Schwartz. The potential is read from 
a 110 point mesh. The number of cores (SCORES) 
and the valence orbitals (WALES) are specified 
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along with the ionicity (ION). The self 
consistent criterion (TOL) and the eigen value 
accuracy criterion (THRESH) were fixed at 0.001 
and 0.00001 in the 441 integration mesh. The 
nuclear charge (Z), the state (KHIZ), the occupancy 
in that particular state (W/ffil) and the trial 
energy for that state (SE) were due to Hermann- 
Skillman [30] , The input data was in a 110 
point mesh. The values at the points not included 
in the presentation mesh are obtained by intei^ 
polation. 

(ii) Generation of molecular potential. 

The program consists of subroutines for Lagrange inter- 
polation and for calculations of co-efficients of IlbVdin's 
alpha expansion [44-46], used to express potential about 
any center and a subroutine for integration by Simpson's 
rule . 

(iii) Non-SC? part 
Input parameter 

(a) The radii of the Watson sphere and the atomic 
spheres 

(b) Charge on the Watson sphere 

(c) Molecular potential (output of (ii)) 

(d) Symmetry inputs This consists of a central 

atom and combinations of ligand atom basis functions, 
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transforming according to a particular irre- 
ducible representation of the cluster symmetry 
group. The angular momentum quantum number and 
the magnetic quantum numiber values along with 
parity are also specified, 

(e) Trial energy value as The deep levels are 

treated as frozen coi’es. The valence orbital 
of the cluster are found by scanning the energy 
in steps of 0.01 Rydbergs. 

The program consists of the following subroutines; 


IKPUT ; Reads and checks the consistency of the input 
data. 

SETUP; Roads and chocks the consistency of the symmetry 
input. Determines the size of the secular miatrioc 
in its' reduced form, using symmetry. 

TMAT and GMAT: Calculates t matrix end G matrix. Radial 

equation is solved for each energy value. Starting 
functions are calculated in subroutine PSTART. 


Sli(TX; Using above two subroutines, completf secular 
matrix is set up, 

EIGEN; Solves the secular matrix and finds eigen values 
and eigen functions. 

OSBP; Calculates ordina,ry and modified spherical Bessel 

function, Hankel function, Neumann function and their 

“ f 

derivatives. 



CGC: Evaluates the Glebs ch~Gordon co-efficients. 

YLM1 ? Computes radia.l spherical harmonics. 

Other subroutines such as IKTSGR, IKTSEP and LlilEQ are used 
for integration, interpolation and solving simultaneous 
linear equations respectively. 

(iv) SCE PsjTt 

In addition to the non-self consistent input the 
occupancy of ea.ch state with the non-self consistent 
energies are specified here. The maocimum number of 
iterations is also specified in case the self consistency 
with an accuracy of 0.001 Ry. is not obtained. In addition 
to the non-self consistent subroutines we have 

MSSCE; Controls iterations and stops program a.fter either 
self consistency is achieved or a specified 
number of iterations are done. 

I'lEliLIZ; Normali2ES the wave functions at every iteration. 

Calculates charge distribution in different region. 

VGSNs Generates potentials for new iterations from 

normalized charge densities of previous iteration. 
It also calculates various contributions to the 
total energy of the cluster. 

The self consistent energy levels and noriiialized 
wavef unct ions thus derived form the basis of the analysis 
made in the present work. 
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CHAPTER II 


^JA^dn-oruc, cAoicdc 

II. 1 Introduction 

Physicists and physical chemists have a sustained 
interest in developing quantitative description of chemical 
bonding of molecules and solids. The systems which have 
been substantially focussed on are the small molecules [1] 
and transition metal complexes [2], Their electronic structures 
nave been determined and the question if their bondings with 
tneir ligands are ionic or covalent have been given careful 
consideration. The present work is an effort to make ab-initio 
tneoretical analysis of the bondings in ionic and covalent 
solids , 

Historically the crystal field model has leen exten- 
sively used to interpret the absorption spectra [3 ] » The 
model proved to be successful as a semi-empirical attempt to 
calculate the crystal field splittings from, pure electro- 
statics [4 ]. However, the experimental results indicate that 
the crystal field approximation is quantitatively very poor, 
being incorrect even to the first order. The active 
interaction of the ligands is brought in 

through the molecular orbital approach. This improved the 
agreements. But, the coefficients of central atom or ligand 
wave functions had to be determined by neutron scattering 
or nuclear magnetic resonance experiments. If the Hartree-Bock 
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equations are solved exactly along with the relativistic 
effects then the theoretical results without parametrization 
are expected to agree sufficiently closely with the experi- 
HBntal results. However, this calculation would take hours 
even with large computers and hence become very expensive. 

We have used the multiple scattering Xq; self consistent method 
which is an approximation to the Hartree-Fock method. 

II. 2 Chemical Bond 

When two or more atoms are brought close together, 
the nuclei and the electrons of these atoms undergo coulomb 
interactions with each other resulting in fixed inter- 
nuclear distances atO°X and a redistribution of the electronic 
charge density, corresponding to the minimum of the total 
energy of the system. The atoms do not maintain their exact 
original identity and, the totality of the nuclei and the 
electrons form a single system. The system is called a 
molecule if a small numbers of atoms v/ere originally involved, 
a complex if it has a total electric charge and a solid if 
the number of original atoms is very large ('^10 ) and the 

equilibrium positions of the noclei form a periodic lattice. 
However, some semblance of the original atoms are maintained 
in that the core electrons around each nucleus usually do 
not change much and what were atoms are now identified as ions 
in different degree of ionization. In equilibrium, each ion 
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is located in a position dictated by the total potential 
miniiEum, A very small perturbation, leading to a harmonic 
approximation, will generate collective oscillations of 
these ions or ps'udoatoms. Hence we may form the classical 
picture of the total system as a collection of ps'udoatoms 
bound with each other by stiff springs. An imaginary spring 
between two neighbouring psludoatoms is called a chemical 
bond . 

As we see, the real mathematical problem is the 
solution of Schrbdinger equation for a system of m electrons 
and n nuclei. Hov/ever, tnis is insoluble when m and n are 
large. Hence different mathematical and physical models 
are made to describe the binding patterns of various systems. 
Of these, three nxDdels are very prominent ; (i) ionic, 

(ii) covalent and (iii) metallic. In each model we have to 
define the psudoatoms, the distribution of the left over 
electrons which is in some degree of delocalization with 
respect to the charge distribution for an atomic orbital, and 
a description of the imaginary entity called the ch^^mical 
bond in terms of the charge density distribution of the 
delocalized electrons. In the present chapter we shall be 
interested in the systems, which are best described either by 
the covalent model, or by a mixture of covalent and ionic 
model. Hence we begin with, a brief description of ionic and 
covalent bonding models. Our description is some what 
coloured by the historical development of these ideas. 
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II. 3 Ionic Bonding 

If we have two kinds of atoms A end B with their 
electronegativity differing greatly, the transfer of an 
electron during their interaction converts them into 
oppositely charged ions 

A = A"^ + e + I^; 

B^e = B Eg 

where I^ and Eg are ionization energy of atom: A and electron 
affinity of atom B, respectively. The electrostatic 
attraction set up between A and B cause the formation of a 
molecules 

A"^ + b" = A^B“ + E 

where E is the energy change. As a result the ions in the 
molecules must be at such a distance that the attraction 
between the electron shells and nuclei is balanced by the 
repulsion between electron shell themselves or between the 
nuclei. Eor a solid, the attraction between the cations and 
the anions is, at the equilibrium distance, balanced by 
the repulsion between the cations or anions themselves. 

We may visualize the cations and anions as psudoatoms and 
as if bound by stiff springs and saj that an ionic bond is 
formed. The ions may be considered as spherical balls, 
since in the process of ionization, the electron shells 
of the atoms acquire the characteristics of 
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the noble gases. This view was developed by Kossel in 1916 . 
However, a complete transition of the electron from one atom 
to anotner to form a pair of cation and anion never takes 
place [ 5 ] there is never a cent percent ionic bond. 

II *4 Covalent Bonding 

Consider a molecule consisting of atoms which have 
same electronegativity, e.g. K2. Here the two atoms are 
equal in their power to attract an electron and hence 
ionization does not take place. In I9I6-I9I8 Lewis and 
Langmuir suggested that the chemical bond is formed by a 
pair of electrons held jointly by two atoms. Such a bond is 
termed as covalent or homopolar. The sharing of pairs of 
electrons is conventionally represented by a line or two 
dots between the symbols of the elements. For example, 

H ; H ; N ; ; ; N ; S : : S ; Cl ; Cl 

H - H ; F = K ; S = S ; Cl - Cl 

The diagrams are called Lev/is structures. The ionic bond 

can be described rather well in tsrrrcs of forces between 
the charged ions, v/hereas the covalent bond requires a 
more detailed recognition of the wave nature of the 
electrons. 

In 1927 Heitler and London showed that if a hydrogen 
molecule consists of Cwo hydrogen atoms a and b in their 
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1s ground state, the e le ctronic wave function a(1 ) b(2) 
gives only a slight minimum for the wrong internuclear 

'j 

distance. Wnereas its symmetric component 2 [&(”' ) b(2) + a(2) b(1)J 
shows a deep energy minimum for approximately correct inter- 
nuclear distance R = . Apparently the stabilization of 

the chemical bond depends on the exchange symmetry of the 
electrons. This successful work of Heitter and London opened 
a nev/ field in quantum chemistry and gave rise to the socalled 
valence Bond theory. 

Soon after this Rulliken came up with the idea that 

the valence electrons of atoms may be considered as belonging 

to the molecule as a whole. This led to the idea of molecular 

orbitals. In this I/O-LGAO scheme, the molecular orbitals are 

formed by linear combination of atomic orbitals available for 

the valence electrons. Bor the case of hydrogen molecule the 

1 

two electron wave function is then given by U = 2 [U(^ ) ■. 

+ U(2)] • Such a wave function is symmetric under 

permutation of the two electronic coordinates and shows a 
similar stabilization of the covalent cnemical bond to that 
of the He it ter- Bond on wave function. Both the wave functions 
have to be supplenanted with a spin factor 

corresponding to the antiparallel spins of the electrons in 
the singlet ground state. The MO idea turned out to be 
exceedingly fruitful, particularly in molecular spectroscopy, 
since the molecular orbitals reflected the basic symmetry 
properties of the molecules involved. 
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The most physical description of the ICAO approxi- 
mation is that the electron moves in an orbit which extends 
to the neighbourhood of both n-,clei. In the neighbourhood 
of one nucleus, e.g., A, the MO resembles an atomic orbital 
U . Similarly in the neighbourhood of the other nucleus B, 
the MO resembles Ui^. Since the complete i'-C has characteristics 
separately possessed by and , it is a natural step to 
adopt the method of linear combinatioris 


U = 


0 u + c, u, 

a a b b 


( 2 . 1 ) 


hj s js most conveniently vncitten as 


TJ = K (U + MU. ) 


( 2 . 2 ) 


where U is the normalizing factor. The value of N is given 

by 




5aa + 2 X 


= 1+2 ^ 


( 2 . 3 : 


when U, and U+ are normalized and 

cl D 


Mb = ^ 'M M d d 


(2.4) 


Por most purposes, the normalisation factor can be dropped 
to obtain 


U = U + M Lb 
a b 


(2.5) 
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The constant X rrsasures the polarity of the orbital and 
may have any value ranging +co to -oo , according to the 
nature of the combining atoms. According to the Ritz 
variational method this value of 7 ^ must be choosen to minimize 
t he e ne rgy fun c t ion 

4tj H U d'C rN 

4 = p (2.6) 

4U dr 

We still have to decide which orbitals and Tj^ of the 
two atoms may be combined together. In order to have an 
effective combination between a given and it is 
ne ce s s ary ; 

(a) that the energies of and in their respective 
atoms should be of comparable magnitude; 

(b) that the charge clouds of Ug and should overlap 
one another as much as possible; 

(c; that Ug and U-j^ should have the same symmetry 
relative to the molecular axis ab. 

It is difficult to choose pseudoatoms in the case of 
covalent bonding as well as w'orking out the forces between 
pseudoatoms since the charge distribution around the nuclear 
site^ are not spherically symmetric. Iwn extreme cases 
that could be conjured up are (i) the pseudoatom is a 
sphere of the same size as the original neutral atom, 

(ii) the pseudoatom is the ion consisting of the nucleus 
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and the core shell electrons. The third model, which is 
the case in between these, has to be defined after the charge 
distribution is worked out. We draw the charge distribution 
contour map, identify the outermost contour surrounding each 
nucleus and define this as the boundary of the pseudoatom. 

This third model is more realistic, and it forms the basis of 
talking of bonding in absolute terms, but is rather difficult 
to handle. Model (i) is simple, allows us to study the 
charge distributions relative to the neutral atoms, i.e., 
the charge transfer to reach the molecular or cluster equili- 
brium with respect to the free atom.s and as such permits a 
simple description of bonding. We may draw the contour 
map of the relative charge distribution and thus have a visual 
description of bonding [6] . Or we may determine the charge 
transfer from each type of pseudoatom, and plot them in a 
charge transfer diagram to get a simple quantitative measure of 
the ionic, covalent or mixed ionic- covalent bonding. Details 
of the charge transfer diagram will be discussed when we take 
up the specific covalent or mixed systems later in this 
chapter. 

II. 5 A Critique of Bonding; Theory 

There is no theoretical foundation for the electron 
pair bond of Lewis and Langmuir. It is not known why two 
and not some other number of electrons should be required 
for a single covalent bond, or how they could act as a bond 
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between tbe atoms and at the same time form part of the 
electronic system of the atoms. Secondly, no explanation 
was offered of the angles between the bonds formed by an 
element. The tetrahedral distribution of four, and the 
octahedral arrangement of six had been deduced for a number 
of elements from classical stereochemical studies, but 
these and other symmetrical arrangement of bonds could not 
be predicted theoretically. Again, considering only the 
valence electrons there vras no way of deciding whether 
a molecule or complex ion should be foimulated with ionic or 
covalent bonds. For example, F and Cl have the same number 
of valence electrons, but shows different chemical properties. 
Also any satisfactory theory of interatomic bonds must take 
account of the deeper structure of the atom, which plays a 
part in influencing the chemical individuality of the atom. 

II. 6 Symmetry and Basis Functions 

The symmetry property of the complex helps in reducing 
the complexity to solve the set of linear equations obtained 
by applying equation 2.6 in the previous section for different 
orbitals. If the point group is at hand, we can construct 
the structurally symmetrized wave functions. These wave 
functions are used as input in the self consistent process 
of the calculation. 

The LCAO method may be used to construct orbitals of 


t he form 
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w - f (v) + 2 (2.7) 

where ( ^) is a wave function of the central atom transforming 
in the molecular or crystal point group as the irreducilsle 
representation t and ^ ^ linear combination of the 

ligand wave functions transforming as the same irreducible 
representation v , Under this condition only, any smearing 
out of the electronic orbitals is possible over the whole 
molecular unit, since only then we can have any cross terms 
between the two parts of the total molecular wave function. 

In this chapter, we concentrate on the octahedral complex 

where six ligands sit symmetrically on the three orthogonal 

carta.sian axes. The c -bonds are MO with rotational symmetry 

about the line joining the nuclei. The combinations of 

atomic orbitals leading to the c bond are; s-s, s-p^, 

p^ - p , p„-d^ etc. -bonds are IvD with nodal plane 

containing the line between the nuclei fp -pjP -p,...). 

The MO of 6-bonds posses two nodal planes (d - d ,. — ) [7 

xy yx 

So there can be s, p^, d 2> orbioals in the o -bonds, 

z 

Figure 2.1 shows that the z-axis is always pointed towards 
the central atom. Hence all the o-bonds are described by 
using the z-coordinate and the ''--bonds by using the x-y 
coordinates. In addition to this, all the ligand coordinate 
systems are left handed. This choice has been made because 
the fact that two center integrate usually are evaluated by 


of 



48 


using one right and one left handed coordinate system [ 3 ]. 
finally, the positive ends of the twelve x and y axes have 
been chosen in such a way that only positive signs appear in 
the linear combinations of the ligand orbitals making up 
the tp,^-7i orbitals. This is not necessary but since the 
most important ic-bond is the ^2g convenient to 

use only plus sign. The six - ligands give rise to six such 
bonds; consequently, six linearly independent orbitals can 
be constructed using these wave functions. 


Table 2.1 combines linearly the ligand functions with <7- 
and TC- bond characters for an octahedral symmetry corres- 
ponding to each irreducible representation of the octahedral 
symmetry group. The central ion orbitals indicates the 
functions with same symmetry behaviour as the linear com.bi- 
nations in the sartE row. The functions in the sane row do 
not have the sans irreducible representation, rather they 
transform as the same column of the appropriate representation 
matrices [3], 

In this chapter we proceed to study the energy levels, 
charge distributions and bonding of metallic fluorides and 
oxides. We scart with some ionic solids to see the effective- 
ness of the cluster model. 




TABLE 2.1 



II *7 Bner,<.y Levels and Charge Distributions 


c ■‘S Q , ' 5 ”1 

^ o o 


II*7.1 Ionic Solids 
II.7 . 1 a Sodium Chloride 

NaCl is the arcnetypal ionic solid of IB type 
(Figure 2.2). The crystal lattice parameter at 1 8 and 26 °C 
are 5.6398 and 5.6406 a.u., respectively. A '''3 XaSCF 
calculation is performed for the 7 atom (Fa"*" 6 Cl ) cluster 
to find out the one electron energy levels and the charges 
within each atomic region. 

There is bunching of valence levels in groups as 
seen in the energy level spectra (Figure 2.3). Similar 
trend is observed by Brescansin et al. [9 ] . Ibe energy 
of the last filled level is 0.64 Ey, which corresponds 
to approximately 8.7 eV. The experimental energy gap which 
is the difference between first unoccupied and last filled 
level is found to be 8.97 eV [13]. 

The charge distribution is compared (in Table 2.2) 
with a theoretical calculation of a cluster as large as 
27 atoms [ 9 ] and also with the experimental data [lo]. 

The excellent agreement wit n the experimental data gives 
confidence in our calculations. 




trs 
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Table 2.2 


Chars© 

Distribution in the 

Different Regions 


of the Cluster of 

NaCl 


Region 

( 

' 1 Ka"^ 6 Cl" 

I 

( 

I Expt . [9 ] 

' 1 Na'^SCl lONa"^ 


: (e) 

i (© ) 

00 

0 

I 

Central ion 

9.9271 

9*9 +4?^ 

10.1984 


(10.1452) 



!Te arest 

17.5362 

17 . 6 + 45 '. 

16.9173 

neighbours (n.n 

) (17.1142) 



Next n.n 

- 


10.1 176 

Next next n.n 

- 


16.7002 

Extramolecular 

0.5073 


0.6198 

region 

(0,7470) 



Intersphere 

2.3481 


14.3477 


(4.4226) 




The bracketed values are for the second lattice parameter at 26 °C. 


The information available in the literature is that 
the lattice parameter of KaCl is 5.6398 a.u. [11 ] and that 
the ionic radii of Ka and 01 in their closed snell 
configuration are 1.96 and 3.62 a.u., respectively [12]. hfe 
note that the sum of the ionic radii is less than the lattice 
parameter. Hence some adjustments are possible in the 
values of the individual -atomic radii keeping their sum 
equal to the lattice parameter. But we know that the ionic 
radius would affect our results strongly since the basis of 
our method is matching the logarithmic derivative of the wave 
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functions at the ionic sphere boundary. Hence we have tried with 
the various combinations of ionic radii (Table 2.3) and 
compared the results with experimental values. It is evident 
that the combinations in which one of the ions has its 

radius equal to, the ionic radius give very close agreement. 

7— 

We choose the/ combination to find out the charges in the 
different region of the cluster. The last set of radii in 
Table 2.3 is obtained by minimising the volume of the 
intersphere region. There is a considerable mismatch with 
the experiments in this case. Thus the concept of minimizing 
the volume of the intersphere region (a concept used for 
linear molecules) at the cost of atomic radius fails, The 
radius of the central ion is also likely to vary with the 
environment it experiences. Hence an error in the choice of 
the radius can by no means be removed, but, the ionic 
radii is a good point to start with. 

Table 2.4 shows the percentage ground state charge 

to 

density population of HaCl cor i-es ponding/each of the valenoe 
orbitals. ^os’ *^1 intersphere 

charge, outersphere charge, central ion charge and the ligand 
ion charge respectively. IR signi-ies the irreducible 
representations of the octahedral point group for our case. 

The intersphere cnarge is found to increase as one 
goes up in energy for a particular IR. This shows that as we 
go up in energy, the orbitals take up more and more covalent 
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Table 2.4 


Per cent age Ground State Charge 
Density Population in PaCl 



d - 5 

.6398 

(180G) 


d 

= 5.64056 (26° 

c) 

IR 

^int : 

q 

^os 

1 q 

i ^1 ^ 

^int 1 
1 1 

^os 

1 q 

1 

t 1 

^1 


0.06 

0.00 

99.93 

0.00 

0.06 

0.00 

99.93 

0.00 


0.34 

0.00 

99.61 

0.00 

0.35 

0.00 

99.61 

0.00 


3.54 

0.14 

0.15 

16.03 

3.53 

0.15 

0.15 

16.03 

2^1u 

2.97 

0.12 

0.09 

16.14 

2.97 

0.13 

0.09 

16.14 

1 

y 

o 

2.88 

0.00 

0.00 

16.19 

2.60 

0.22 

0.06 

16.19 


1 5.55 

1 .40 

3.12 

13.32 

15.47 

1 .46 

3.15 

13.32 


19.04 

0.00 

0.00 

13.49 

16.96 

1 .51 

0.10 

13.57 


5,31 

4.54 

0.59 

14.87 

6.41 

4.42 

0.76 

14.73 


13.75 

0.00 

0.00 

14.37 

12.54 

0.09 

0.04 

14.42 

'"lg 

10.88 

0.00 

0.00 

14.85 

10.18 

0.34 

0.00 

14.88 

2e 

g 

12.53 

0.00 

0.00 

14.58 

7.78 

3.43 

0.58 

14.70 


1 1 .07 

1 .93 

0.34 

14.44 

11 .03 

1 .98 

0.34 

14.44 
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character. The deep valence levels, e.g., la. are hardly 
affected by the sli^;ht change in the lattice parameter. 

But lor the shallow valence levels, e.g., 2e^, the charge 
in the ligand increases slightly at the cost of the inter- 
sphere region. This indicates that more of the atomic 
character is maintained at lai'ger lattice parameter, 

II, 7. 1b LiB and Baf 

Both LiF and NaP are knovm to have a structure of 
NaCl type (Figure 2.2). Table 2.5 indicates the numerical 
values of the alkali-metal-ligand distance (d), the Watson 
spuere radius (r ), the central ion radius (r ), ligand 
radius and the exchange parameters (o^), used in our 

MS Xq,SCB calculations. 


Table 2.5 

Numerical Values of Different .Parameters Used in Our 
Calculat ions for Lil' and FaF 



^ $ 

d(a,u. ) 1 
» 

t 

r^(a. u. ) 

I r^(a.u. ) 

t 

r-j (a.u. ) ) 

C5 

I'x 

1 a 
: ^ 

LiP 

4.06 

6 .46 

1 .66 

2.40 

0.73087 

0.7365 

BaP 

4 . 62 

6.58 

2.66 

1 .96 

0.73044 

0.7365 


The cc values used in our calculations are due to Schwartz [14]. 




<^8 


Figure 2.4 shows the energy level diagram of liF. 
The feature tnat imnediately attracts attention is the 
bunching of the valence levels in tyo groups. Evidently 
the system is ionic with a small amount of coavalency . The 
eigen functions of the upper set are fluorine 2p-like and 
those of the lower set are fluorine is-liJce. The small 
dispersion in each set is due to covalency effects, the 
addition with a small coefficient of the Li orbital of 
symmetry appropriate to the fluorine combination described 
by the irreducible representation of the level. NaF energy 
levels also snow similar features. The X-ray photoelectron 
spectrum of LiF finds the F 2p and F 2s energy difference 
to be 1.54 Ry [15]. Our calculated shift is 1.43 Ry. 

The self consistem:ly obtained charges in the 
different region of the cluster are shown in Table 2.6. 


Table 2.6 

Self Consistent Charges in ^ Lifferent Regiongof the 
Clusters of LiF and I’laF m units of e. 


I 

Region i 

(LiF) 5- I 

(Fai! y 1 

(FaF)"^ 

Extramole cular 

5.005 

1 .743 

2.344 

Central ion 

1 .993 

10.226 

9.761 

Ligand ion 

9.399 

8.844 

8.511 

Intersphere 

0.609 

4.968 

5.823 

Watson sphere charge 

5.0 

5.0 

4.0 

Constant Potential ' 

-0.654 

-0.359 

-0.514 


lergy (Ry) 
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be less than the true ionic distribution, shows that the 
ionic part of the mixed bonding may be in reality "super" 
ionic where ionicity is overdone . To look further into 

a theoretical possibility of superionicity , we have soiled 
the MS Xq-SCP for a cluster of 1 Na”'"'*' central ion and 6 S' 
ligand ions column 3 of Table 2.6. see that such a 

system stabilizes by coming closer to the DTa charge distri- 
bution and transferring more charge to the intersphere region 

5— 

compared to the standard (KaP) case. This means that if 
one electron of a Na atom is knocked off by a suitable 
radiation, the local environment will respond by a redis- 
tribution of charges corresponding to more covalent character. 
However, it is to be noted that the charge within the 
ligand ions i-s not much affected. But the energy level 
structure is substantially changed. 

Table 2.7 shows that the a-^ g remain 

core levels of Na"'' for all practical purposes. levels 

o 

and 1e remain approximately P core levels. The level 2t^^ 
shows a strong bonding between the cation and the anion. 

All the higher levels have similar charge distribution 
pattern; almost empty cation region and large ring like 
distribution in the anion-annulus, the density slightly 
stepped up at the boundary as we go from the intersphere 
region into the anion. 
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Table 2.7 


Perce nta^:e 

Ground State Population Density Distribution 

Por (NaPg) 

^ Cluster 



] 

IR 1 

1 

q. 

^os 

1 

! "1 

I 

j '^int 


0.00 

99.94 

0.00 

0.06 


0.01 

91 .39 

7.26 

0.86 


0.21 

0.22 

1 5.80 

4.78 

1e 

g 

0.33 

0.15 

15.88 

4.21 

2^1a 

0.17 

8.38 

14.58 

3.99 

3a^g 

1 .08 

3.42 

13.73 

13.13 

1t^ 

2g 

1 .25 

0.23 

13.95 

14.82 


0.72 

0.05 

14.52 

12.25 

2®g 

2.44 

1.13 

14.73 

8.01 


1 .24 

0.65 

14.58 

10.61 

^*ig 

0,66 

0.02 

14.80 

1 0.71 

'’■''lu 

3.12 

1 .06 

14.87 

6,59 
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II. 7. 2 Transition Metal Oompounds 
11,7® 2a Tluorides of Transition Ivatal 

shall discuss the different transition metal 
compounds like lIFi?^ , FeP 2 j Cokg? Nil '2 ZnF^? 

which are known to have covalent character. The first 
two of these transition metal compounds have a perovskite 
structure (Figure 2,5a) while the others have a rutile 
structure (Figure 2,51). In the case of perovskites the 
central atom is a transition metal with six equivalent 
ligand atoms placed at the face centers of the unit cell. 

The K atoms occupy the cube corners. We do not consider 
the k atoms explicitly. But their stabilizing effect is 
taken care of by putting appropriate charges on the Watson 
sphere. For the IFg (k = Fe^"^, 00^“^, and Zn^*^) 

complexes, we find four of the six F atoms to be equidistant 
from the center while the other two are slightly away. This, 
tnus, distorts the regular octahedral symmetry, generating 
an axial field. Again, the four equidistant atoms are not at 
right angles but are at an angle slightly less than 90°, 

This asymmetry lifts the planer degeneracy in energy. 

The a --values of Fe , Go, Fi, Zn and F were taken 

from Schwartz [U ] to be 0.71151 , 0.71018, 0.70896, 0.70673 

and 0 , 73732 , respectively. The values of the radii used 

(Table 2.8) are obtained from 

r . 

r. = d X - — I i == c,l (2,8) 

cl 




\ X ' ' 


ov 


pi:; ■' t s ?' 1 f K i'/f F ■:» 



FIG. 2.5 b 

sfruciura of a rutfle MF 2 



65 


\.;here i is the type of ion considered , c and 1 denote the 
central and the ligand iore respectively and d is the metal- 

ligand distance. 


Table 2.8 


The Value oi Radii Used lor Different TM Compounds 
in Atomic Units 


Compounds 

5 

1 

1 

1 

! 

I 

^0 

« I 

5 i 

I 

^1 

i 

i ! 

1 r^ 

1 1 

i 

la^eU^ 


5»617 

2.653 

1 .432 

- 

U-eip 


5.258 

2.525 

1 .365 

- 

itr-iif^ 


5.129 

2.456 

1 .336 




5.396 

2 .584 

1 .406 

- 

ppP 

ie- 2 


5.675 

2,729 

1 .331 

1 ,473 

CM 

O 

o 


5.560 

2.640 

1 .460 

1 .436 

Uil2 


5.571 

2.570 

1 .510 

1 ,389 

Z 


5.536 

2.624 

1.456 

1 .429 

The two values 

of 

KPel’.v and 

3 

KFiP-, correspond to the 

3 

eases of 

different central 

ion-ligand 

ion distance 



The atomic 

ten-fold 

degenerate 3d 

-level breaks up 

into a six-fold 

degenerate t 

„ level and 
2g 

a four-fold 

de gene rat 


e^-level in presence of the crystal field. The e^- level has a 
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lower energy as compared to the t 2 ^-level in an octahedral 
lie Id. The difference in energy between these txv'o levels, 
called the crystal field splitting, is designated by 10 Dq. 
Figure 2.6 presents the one electron spectra of FeF 2 ? CoFg 
and F 1 F 2 •> Tbs levels are seen to split due to the asymmetry 
in their structure. The XfS experiment [15] findsthe P 2s 
and P 2p energy difference for FeP 2 to be 1 ,54 Ry which is 
in very good agreement with our calculated value of 1.36 
Ry, The small discrepancy is due to the correlation effects. 
However, before we start estimating correlation energies of 
these levels, we should have experimental data on more 
transitions as also got our calculations more refined to 
take into account at least the effects of the second nearest 
neighbours. The 3d occupancy and the crystal field 
splittings in the different transition metal compounds are 
tabulated below (Table 2,9). Our results for the perovskites 
compared withtJuase of Adachi et al. [17 ] , 

The rather good over all agreement for the 
perovskites shows that the splittings of the d -leva Is into 
t 2 o. and e^ should not be thought of as due to the influence 
of a perturbing electric field as is done in the simple 
crystal field theory. The only resemblance that the 
muifir-tin potential used in our approach has the electro- 
static potential of the crystal field theory is the 
octahedral symmetry. The splitting in our case depends 
crucially on the fact that we are using molecular orbitals 


are 



Relative energy (Ry) 



FeF 2 C0F2 NIF2 


Fig. ' 2-6 


. energy In Ry for FeF2*CoF2 qnd NIF2 • 
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Table 2 ,9 

3d Occupancy and Crystal Held Splittings in the Different 
Transition I'fetal Compounds 


Compounds j 

Present work (e) 

Adachi 

10 Pq 

I Experiment 

1 


et al.[l7](e) 

(present 
work eV) 

(eV) 

KTePj 

5.287 

6.147 

0.83 


KlTiP^ 

3 

7.715 

8.199 

0.88 

0 . 9 [ 18 ] 

PePp 

5.331 

- 

0.95 


GoPp 

6.389 


0.68 


MPp 

7.956 

- 

0^61 


ZnP2 

9.901 


— 



which are solutions of the same effective Hamiltonian and 
occupied in accordance with the Pauli exclusion principle [19]. 

The charge distributions in the regions of the 
cluster calculated self consistem ly are shown in 

Table 2*10« The lattice parameters are 4.085 and 3.893 a.u. 
respectively for KPeP^ while it ia 3.793 and 3.99© a.u. 
respectively for KlJiP^. 
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Table 2.10 


Ctiarge Distribution in the Different Regions of the Cluster 
for the Transition Metal Compounds in units of e . 


Complexes 

i i 

1 Extramole Gular j 
1 1 

I 

[ Gantral-ion ; 

I t 

ligand-ion 

1 

! Intersphere 

( 

I 


3o119 

23.940 

7.741 

10.494 


2.574 

23.927 

7.343 

13.440 

3 

3.031 

24.486 

7.547 

13.201 

KNiF^ 

3.022 

26.816 

7.593 

10.603 

Pep2 

3.700 

25.01s 

7.141 

12.438 

C\) 

o 

o 

3.105 

26.029 

7.503 

10.847 

PiP2 

0.599 

26.262 

8 . HI 

10.299 

ZnP2 

1 .008 

29.622 

7.940 • 

10.076 


The expected ■general trend in both the structures is 
a gradual increase of the total central ion charge as Z 
increases, in conformity it h the simple chemical picture 
based on the Pauli electronegativity scale. Also, as the 
interatomic distance increases, there is less mixing of 
the ligand electi-ons with that of the central-ion and the 
charges within the central- ion and the ligand-ion region 
are expected to be atomic like. The casesof Ki^eP^ and 
KGoP^ (see Table 2.10) clearly point this out. The charges 
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in both the regions increase with distance creating more 
of atomic environment around tfat^^ir ovm centers* 

The basis functions used in our calculation are 
the same as those used as molecular orbitals in the ligand 
field theory* Hence their coefficients, derived from our 
self consistent calculations, are the same as the bonding 
parameters \^and/or b-j-r * In the molecular orbital theory, 
and /v^ can not be calculated from first principles. They 
have to be deduced by matching the parametrized orbitals 
form factors with neutron scattering experimental data or 
their spin-density at the nucleus with NFE experim.ental data. 
There are not many experimental data available. But neutron 
diffraction data for KNiB^ and KBeB^ do exist [ 20 ], and they 
are 5.11 percent and 6.24 percent^ respe ct ive ly. Our calcula- 
tions give ^0 efual to 7.55 percent and 5.58 percent for 
KNiB^ and KBeB^ ^respect ively. The agreement is quite satis- 
factory. lAfe may note that the neutron scattering data are 
due to orbitals with unpaired spin only and hence a spin 
unrestricted calculation is expected to give still better 
agreements. 
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II»7»2b Oxides of Transition Fietals 

fO'onoxides of the 3d transition elements form a 
very interesting group of compounds with a wide spectrum of 
electrical and magnetic properties. In this comparatively 
small groupj ferromagnetics and antiferromagnetics^ dielectrics 
and compounds with high metallic conductivity can be found. 

Such an 'unusual combination of properties in a series of 
iso- struct Ural compounds suggests very strong differences in 
electronic structures, chemical bonding and energy spectra of 
compounds. Therefore, it is not surprising that for a 
long time the electronic structures of the 3d monoxides 
has been the subject of a large number of experimental and 
theoretical studies [ 19 , 21 ], 

II.7„2b.1 MgO 

The one electron energy spectra of MgO is shown 
in Figure 2,7» ^ see ntliat the set of occupied 

orbitals bet'v.een -0,56 By and -0.68 By have a predominen- 
tly ligand 2p character. The energy of the last filled 
level is 7.62 which is in very good agreement with the 
experimental energy gap of 7.8 eV [22], The bonding 
effects between 1% and 0 ions and among 0 ions split the 
valence band into tv/o groups. One group of 
It. , It,, , Ito.., 3 t 9 „ and 2a-, levels is the set of the 
non-bonding type levels with nearly complete 0 2p character. 
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Tne other group includes 1 s 2t> ^ 1e^ levels^ The wave 

1 g I u. g 

function of the la., , 2t. and. 1e levels correspond to 

I g I a g 

the ligand 2s orbitals mixed with the Mg 3s orbitals» A 
DY T (X calculation was performed by Satoko et al« [23] » 

The energy level structure shows a similar seq.uencfi except 
for the inversion of 3t.|^ and self- 

consistent orbital energies are also compared with those of 
Satoko et al, [23] (Table 2J1), Any difference encountered 
could be due to the different lattice parameters used. Our 
lattice constant is 4.2 a.u. in contrast to 3.98 a.u. of 
Satoko et al. The IPS eoniirins the Mg 2s and Mg 2p levels ['^ 5 ] 
at 6.69 and 3 <,84 Ry, respectively. They axe in fair agi-eement 
with our calculated values of 5.8 and 3.42 Ry^ respectively. 

Table 2,11 

Self Consistent Orbital Bnerries in Ry for MIgO 


Ions 

5 ^ 

5 Orbital J 

T ! 

1 

1 

Present work i 

t 

t 

Satoko et al.[ 23 

Mg 

1 s 

92,2 

CO 

— 

cn 

0 

1 s 

37,37 

37.13 

% 

2s 

5.80 

5.37 

Mg 

2p 

3.42 

3,02 

0 

2s 

1 .67 

1 .49 

0 

2p 

0,51 

0,51 
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If the MgO yere a perfect ionic crystal, the charges 
on the Mg and 0 ions would equal +2e and -ie, respectively. 
Table 2.12 shows the self consistently obtained charge distri- 
butions in the different regions of the cluster. Because of 
the partial covalency, however, sohb amount of electrons flow 
back into the cation region from the anion region. The net 
charge transferred from the ligand ion is 1 ,86e, while that 
from the central ion is 1 .65e. One can assume that the central 
ion looses 2e and then an amount of charge equivalent 0.35© 
flow back. Since the anions only give away electrons, the 
charge flowing back into the cation region can be thought of 
from the ligands. Thus MgO has an inverse ionic character 
(discussed in detail in Section 11,8). 

TjffllE 2.12 

T" 


Charge Distribution in the Different Regions of the Cluster 
of MgO in units of e 


Region 

1 

o 

0 

1 

Central-ion region 

10.35 

10,33 

ligand- ion 

6,14 

6 » 06 

Intersphere 

18,36 


Extramole cular 

4.47 

~ 

Net charge transfer 

from Mg 1 ,65 

1 .67 

Net charge transfer 

from 0 1 .86 

1 .94 
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The covalency parameter, as obtained from 

[24], is has an excellent agreement with 

our calculated value of 6.6755, 

The DV Xq; method does not approxiciate the potential 
within each atomic region to be in the muffin-tin form as 
in hS Xq. method. But apain, to obtain a better result with 
DV method one needs a large number of sample points. To 
conclude, we say tha.t even with this hS X^, method one 
can obtain sufficiently good results so long as the potential 
is not highly asymmetric, 

II, 7,2b. 2 VO and FfaO 

Vanadium monoxide and Manganese monoxide have a 
simple crystal lattice of Fa Cl type, in which each metal 
atom, in the defect free lattice, is surrounded by a regular 
octahedron of oxygen atoms. The length of the retal-oxygen 
bond are set at 4.48 and 4.20 a.u., respectively [11]. 

In accordance with the formal valence of the elements and 
the filling of the energy levels of the compound, the 
octahedral cluster in VO and InO should contain 10 excess 
electrons as in the case of MgO. Tke V and Mn radius were 
set at their muffin-tin value of 2.40 and 2.31 a.u. [ 25 ], 
respectively. The c;_yalues of V, Im and 0 were sat at 
0.71556, 0.71279 and 0.74447, re spe ct ive ly [ 1 4 ]. 
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ri ^ 1 

The ener9:y level diagrams corresponding to (7'Og) 

1 0 — 

and (MnOg) clusters are shown in Figure 2^8, The results 

1 0 “* 

are plotted in a relative scale » shows a similar 

clustering of levels as in the case of !g0 whereas in hnO 
the upper levels are somewhat spread out. The 0 2p level lies 
6.12 eV below the tg^^. level (Fermi level) which is in contrast 
with the available XF3 data with an energy 4 eV [26]. for MnO. 


Table 2,13 show's the charge distribution of VO and 

ivihO. Our calculated values for VO are compared witn that of 

Gubanov et al, [27,28 j for a different nstal- ligand distance 

1 - 

of 4.08 a.u. The table also includes the (hnO^) cluster 
calculation results of' Johnson [29] which has a tetrahedral 
syraiiietry. 


Table 2.13 

Charge Distribution in VO and IfoO Clusters (in units of e) 


Region 1 

1 

(V06)'°" 

' / VO 1 

1 iVOg) 

; [27] 

[ (VO. 

: [28: ; 

1 (itoOg)^® 

1 

|(bhO^) 

: [ 29 ] 

Central-ion 

22.121 

21 .884 

21 .91 7 

24,362 

22,462 

Ligand “ion 

7.023 

6,805 

6.809 

7.370 

6.274 

Int ensphere 

7 .262 

13.573 

T)3.689 

11.248 

9.456 

Extramole cular 

9.423 

4.721 

4.535 

ro 

ro 

0 

0.968 


Rflotive energy feV) 
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Fig 2-8 

Om etfctron energy ^ctro of MnO ond VO if^ eV in rcSolive scale. 
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The present work shows a loss of charge equal to 

0.879e and 0»638e from V and fli as compared to the value 

of +2e as considered to start with. Alternatively, we can 

think that the ligand charge has contributed partly to the 

central -ion region and a ma;]or part to the intersphere and 

1 - 

extra^iole cular regions. In cluster both ifin and 0 

loose charge to the intersphere region. For a purely ionic 
case the net charge on V or Fn should be +2e with -2e charge 
on 0 ion. So a cnarge of l.llle and 1 .362e arrive at the 

2 -r 2-h 

ion a.nd I'ki ion. The results of Gubanov et al. shows 

an increase of 0.884 only for the (VOg) cluster. The 

2 - 

charge deviation from 0 ion in our case are 2.97e 

and 2.63e for VO and MnO while that of Gubanov et al. is 

3.19© for VO and that of Johnson is 3.72e. The matching 
”1 0 ^ 

of (FinOg) and (FtnO^) cluster results are never expected 
since the symmetry of the later is tetrahedral. This gives 
a feeling about the structure dependence of the charge 
distribution. In (MnO^)"' cluster, the Fn is in Fin"^”^ state. 

So a charge of 22.462 - (25-7) = 4.462e is pulled into the 
Fn region out of ( 10-6. 276 )x4 = I4.896e of the ligands. 

Since (FinOg) cluster has six nearest neighbours, the 

2 - 

total charge given out by all the 0 ions is (2.63x6=) I5.78e 

of which only 1 ,36e goes to the central ion. This shows that 

(FIqO^) ” is more inverse ionic (discussed in II. 8) as 

1 0 ”“ 

compared to (MnOg) 
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II^7»2b,3 BaTiO^ and SrTiO„ 

l-Z. 2 

We repeated the calculations for the electronic 
structures of BaTiO^ and SrTiO^ for utilization in the 
charge transfer diagram later in connection with the discu- 
ssion of the bonding characteristics of covalent solids » 

Our results agree closely with the existing calculations 
[30,31 ]. 


II»7®3 .Alloys ; Pe-Te Cluster 

The present method may be easily extended to the 
vast area of the localized electronic structure in alloys* 

The alloys usually have short-range structural order. le may 
form a cluster of such a locally ordered unit and work out 
the MS Xq,- SCI energy states and charge distribution. These 
may be used for the determination of a physical property of 
localized nature, e.g. the electric field gradient at the 
nucleus of a Mdssbauer atom. Here we initiate the calculation 
for le-Fi-Te ternary alloy as an example. 

Experimental MSssbauer studies of Pe-Fi-Te alloys 
have been done in detail in our laboratory [32], The 
crystal structure changes at several critical concentrations 
of the netal-metalloid ratio. The typical structures are 
tetragonal, orthorhombic, rhombohedral and defect Fi-As, 

The spectra could be analyzed in terms of two or three or 
four puadrupole doublets with isomer shifts. Each doublet 
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was linked i^ith an le-site in the crystal and the QS and 
IS data were analyzed to obtain detailed infornation about 
the fii'stj second and the third neighbour environnent of 
this le atom, and the complex bonding in this cluster. The 
most dominant bonding was the Pe-Te covalent bonding. Each 
Ee atom is octahedrally coordinated with six Te atoms and 
as such forms a cluster with octahedral symmetry as studied 
in the present chapter, 

for Ee iTe ^ alloy (Ei = 0) the Ee-Te distance is 
5,24 atomic units [33], The Ee and Te radii are 2,23 and 
3,01 a.u. respectively with the Watson sphere radius 8.2 5 
a.u. As initial inputs we put neutral Ee and Te atoms with 
exchange parameters OsYllSl and 0,7000 respectively. The 
electrical neutrality of the cluster is maintained by putting 
zero charge on the Watson sphere. The calculated charge 
distribution in the cluster is tabulated below (Table 2.14), 


Table 2.14 

Cfaar 3 :e Distribution in Ee-Te Gluster 
in unit s of e . 


Sxtramolecular region 

1 .3586 

Central Ee region 

23.8292 

ligand Te region 

50.2941 

Intersphere region 

11 .0144 
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The charge in the central region indicate that Fe 

2 + 

pseudoatom has approKimately the Fe character. The 
li/To'ssbauer study in our laboratory [3 2] shoxv that Fe atoms 
form layers in the crystal and confirm that iron as in 
Fe^"'' state. 

The occupied one electron energy level structure 
of Fe-Te cluster is shown in Figure 2.9. It is seen that 
the intermediate levels have a predominantly Te 5p 
cnaracter . 

The charge transfer diagram (Figure 2.11) to be 
discussed later, show that Fe-Te bonding is strongly 
covalent . 

II , 3 IVfeasure of Go valency 

For charge transfer analysis v/e divide the molecule 
or the cluster in a solid into several regions: One atom 
representing each type of atom or ion in the molecule or 
cluster, e .g . , hydrogenic pseudoatom and fluorine 
pseudoatom in HF, and an interatomic region. The basis of 
the charge transfer analysis is the fact that the neutral 
atoms, when brought close to each other, reorganize their 
intern uolear distance and redistribute their charges to 
come to energy minimum and the corie spending electrostatic 
1 -i nm . In our cluster calculation we have used the 
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(N -Tf) 

Fig. 2 % 

Ont clfctron tntrgy Itvil of Ft-T# 

clufitr. 
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inter-nuclear distances obtained experiL'entally, that is 
the equilibrium distances. Hence our calculation determines 
the equilibrium charge distribution in the different regions. 
Contrasting xs'ith the neutral atom charges in each region 
v/e can find how much charge is flovm in or out of each region 
to reach electrostatic equilibrium. If the charge gone out 
of the cation is exactly equal to the charge entering the 
anion, the system is obviously ionic, iiny deviation from 
this pro csss involves charge being transferred to the inter- 
atomic space, a charge belonging to both the ions or pseudo- 
atoms and hence covalent in nature. Hence from the charge 
transfer diagram we may get a measure of the covalent bonding. 
Probably a more satisfactory measure of covalent bonding 
x/ould be to (use model (iii) as discussed in II. 4) calculate 
the total electrostatic force on each pseudoatom for various 
interatomic distances and calculate the constant cor the 
imaginary spring binding the pseudoatoms. But then the 
calculations would be far more elaborate than the present 
one. Vfe use the charge transfer diagram to analyze our 
d at a ® 

To construct the charge transfer diagram, we note 
that the available electrons in the cluster are solely 
from the central ion and the ligand ions. Since the total 
number of electrons is conserved, we can write 

APc + + ^^ei 


0 


(2,9) 
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I'/here .v P ^ Is the charge transferred to or from the central 
pseudoatom in forming the equilibrium cluster from the 
neutral atoms^ ^ charge transferred to or from 

the ligand pseudoatom in forming the cluster from the neutral 
atoms, APp>i i® charge in the extramole cular and 

int ensphere region in equilibrium state of the cluster. 

Equation (2J) is the equation of a plane passing 
through the origin. The general equation of a plane is 

ax + by + cz + d = 0 ( i , 1 0 ) 

If d = 0 and a = b - c = 1j then the plane passes through 
the origin making an angle 45° vjith all the three coordinate 
axes. The angle, thus depends on the ratios of the 
coefficients. In our charge transfer diagram this depends on 
the charge transfer ratio between the cations and the anions. 
The charge transfer ratio is again a function of the electro- 
positivity and electronegativity of the anions and cations 
and also the structure of the complex. 

Since only two of the three variables and 

yqp . are independent, the description of the charge transfer 
process could be considered complete with the knowledge of any 
two of them® Ife plot our data in the aP^ and APq space. 

The intersection of the charge transfer plane in three 
dine ns ions with the AP^ " APq_ pla^:® is the straight line 
dealing with the charge transfer between the anion and the 
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cation only and correspond to the ionic case of various 
kinds. The amount of charge transferred from the cation to 
the anions is governed by their electropositivity and 
electronegativity, respectively. The line from the origin 
to the ideal ionic point = +1; = -1 for NaCl), 

is named as the true ionic region (figures 2.10 and 2.11). 

The extension of this line in the same quadrant is termed 
as super ionic region. If the electrons are taken up by 
the cations from the anions then we call it an inverse ionic 
system. The inverse ionic line is an extension of the true 
ionic line in the opposite quadrant. The inverse ionic 
character is seen in III-V semiconductors where the ligands 
first give away one electron to the central ion and redis- 
tribution takes place for reaching the equilibrium, for 
the solids considered by us, there are two types of ionic 
lines, one for the charge transfer ratio 1:1 (e.g., FaCl, 

¥igO) and the other for the charge transfer ratio 2:1 (e.g., 

f ef 2 , Wif 2 ) » 

Most systems are mixed ionic and covalent and as 
such their charge transfer points fall off the ionic line. 

If any point is joined to the center of the coordinate system 
then the length v" x) j called the root laean 

squared charge transfer, gives a neasure of the total charge 
transfer for the charge redistribution to take place. Again 
we may determine the coordinates of this point along the 
ionic line and normal to the ionic line. The foimer gives the 



86 


root mean charge transfer between the cation and anion only 
and may be called the ionic bond strength wrhereas the latter 
is the root mean square charge transfer to the interatomic 
and extramole cular region from the central and the ligand 
atoms and is a measure of covalent bond strength. The 
unit of the each bond strength is an electron. Both ionic 
and cova.lent bond strength can also be described as a 
combination of charge transfer ( l^P l\Pj_) - 

Figure 2.10 shows that NaCl has the shortest line 
normal to the ionic line, identifying its basic ionic and 
marginally covalent character. LiF and NaF sho^i? similar 
characteristics, with somewhat more augmented covalencies. 

The transition metal perovskites and rutiles could be 
considered as typical covalent compounds. The ionic bond of 
each of them is much smaller than its ideal ionic limit. But 
the covalent bond strength of each is distinctly large, 
being in each case substantially greater than its ionic bond 
strength. For the rutiles we note that the covalent bond 
strength decreases as we go from Fe to Zn. It is also 
evident that the ratio of the covalent to ionic bond strength 
de crea.ses as we go from Fe to Zn. Hence a,ccording to the 
present analysis, FeF 2 is the archetype of a covalent 
compound with large covalent bond strength in absolute terms 
as well as in relation to the ionicity. The ionic part of 
ZnF 2 shows an inverse character, though small in value. 
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Table 2,15 shovs the numerical values of the covalent bond 
strength of all the systems studied in the present work. 

Table 2,15 

Covalent Bond Strength for Different Systems 
in Units of Elect r ons 


NaCl 

0.2 

KPePj 

o 

C\J 

PePg 

liP 

0.3 

KBiP^ 

00 

CoPg 

NaP 

0.65 



KiP2 





ZnPo 


2.09 

I'lgQ 

2 .51 

Si 

2 .40 

1 .78 

VO 

1.33 

InP 

3.18 

1 .52 

knO 

_0 .85 

GaAs 

1 .13 

1 .12 

BaTiO^ 

3.25 




SrTiO^ 

3 

2.83 




In Figure 2.11 the results of our calculations for 
the oxides of metals aare plotted. I«fe also include here 
the data for a few semiconductors worked out in the next 
chapter^ All the systems show very large absolute values of 
the covalent bond strength. HfeiOj BaTiO^ and SrTiO^ show 
^normal' ionic bond strengths of different degrees. VO is, 
again, a truely covalent system with almost zero ionic bond 
strength. f/lgO, which is known to have semiconducting 
properties, lies in the inverse ionic region. Other semi- 
conductors like Si, InP [34] and GaAs [35 ] also lie in the 
inverse ionic region. The case of IrJ and GaAs were expected to 
be so, since they start vlth a charge transfer from the 
group V to group III elements. Bor Si this is partly due to 
the fact that the central atom and the ligand atom in the 




F'^2.10 

Chpt^. troHslw d toy o m for ^Cl ondf ftuorxks elf metols. 




oroe tronsfer dioqram for cwides of metals, semiconductors an 
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cluster are treated as different, though in real solid no 
two Si atoms could he considered as different in their 
properties. However, the fact that energes is that, within 
the framework of the present analysis, the semiconductors 
could be categorized as a group of substances v/ith large 
covalency combined with Inverse ionicity* 

The present work makes an effort for the analysis of 
the structure of a large number of solids, showing that most 
systems are mixed ionic and covalent, with KaCl in one 
extreme with large ionic and marginal covalent character and 
the semiconductors at the other extreme, with large covalent 
and inverse ionic character. This analysis bridges the gap 
between the purely ionic theories and the band structure 
calculations . 

II„9 Isomer Shift 

Isomer shift is the manifestation of a very specific 
feature of the hyperfine interaction between nuclear and 
electronic charges, namely the finite nuclear size on 
Coulomb energy® The interaction between total electronic 
charge ep„(r) and the electrostatic potential generated by 

nucleons, is 

= -e / pQ (r) \(i') ^ 


U 


( 2.11 ) 
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Por a point nucleus ™ (2.12) 

2 -1 
and U = - Ze / p_(r) V (r) d'c = -Ze <r > 

P © n 6 

( 2 . 13 ) 

-1 

where <r is the average over all electrons present. 

Since U depends only on Z, the energy of a,ll nuclear states 
of a certain nuclide will be same. But for a finite nuclear 
size, the Coulomb interaction between the electronic charge 
and the nuclear potential will be different for the ground 
state and for the excited state. This is indicated by the 
observed energy shift in the Ifdssbauer spectrum. Considering 
the nuclear charge to be uniformly distributed over the 
nuclear sphere of radius R we have 


V^(r) 


& 

E 


[i 


2 J 


r < R 


( 2 . 14 ) 


2e 

r 


r > E 


( 2 . 15 ) 


So the change in potential from that of a point nucleus is 


V (r ) 
n^ ' 


Ze r 3 „ 1 /r%2 

L p 2 ^R'^ 


R 


R 


] r < R 


( 2 . 16 ) 


= 0 


r > R 


( 2 . 17 ) 


The energy shift is now given by 


A U 


-e/ Pq (r) A v^(r) d^: 


( 2 . 18 ) 
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Assuming Pg(r) to be constant over the nuclear volume, 'we 
obtain by integrating the above expression 

A V = I Tt Ze^ P (0) (2.19) 

n 5 e 

Such differences v/ere encountered long ago by the experimen- 
talists in optica,! spectroscopy as isotope shift. To obtain 
the isotope shift one measures the optical transition for two 
different isotopes of an element. The Coulomb energy, which 
is part of the transition energy, v/ill change slightly for 
the two isotopes since R changes with the change in nuclear 
number. 

Now if the source and the absorber consists of two 

chemically different materials then Pg(0) will be different 

for the source and absorber. Hence the fina.1 expression 

absorber 

for isomer shift is then given by (nP^fO)= Pq(0) 

_ p^^Qpource^ 

IS = I Ze^ (-R)^nP^(O) =a Pp(0)- (2.20) 

^ 5 Ey e e 

where ( A R)^ = (R^ ^^^fsrs to two isomeric states 

of a given nucleus. If the relent ivistic corrections a.re 

made the above expression for IS is then multiplied by 

57 

S'(Z) which is 1.29 for Re. 

In our transition nBtal compounds like EReR^, ReRg 
and Re-Te alloy, the isoner shift is calculated theoretically 
using the above expression. Since in all these experim.ents 
^'^Re is used as the source Pg(0) of metallic iron is needed. 
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This value was taken from Walker et al.[36]. The calibration 
constant ct has the term (AR) was choosen to be 0»5'1 [3 6] 
and 0«12 [37], respectively. Table 2.16 shows the charge 
density at the origin and the calculated isonBr shift as 
compared to the experimental data® 


Table 2J6 


Charge Density at the Origin Used to Calculate 
Isomer Shift and Compared with Experimental Data 


System 

1 

p (0) (a.u. ) 

1 

‘ a 

1 ( a, u. mm/sec) 

1 

r— — — — ^ 

t ? 

1 IS ( mm/ se c ) -5 

1 » 

\ 1 

S' (Z)IS 

1 

] Expt , 

1 (mm/seo) 

3 

6,5232 

0,51 

1 .008 

1 ,301 

1 .38 [ 33 ] 

PeR2 

6.5971 

0,51 

0.971 

1 ,253 

1 . 36 [ 38 ] 

Pe “Te 

3.9670 

0,12 

0,544 

0.701 

0.518 p2] 


It may be noted that the value of a as tabulated in the 
literature ranges from 0.11 to 0.51 a.u. mm/sec. The only 
way of finding a better a from experiments is muonic X-ray 
[39] which is drawing attention these days. The calculated 
results show satisfactory agreement with experiments. Our 
results could be improved by taking into consideration the 
effect of the second nearest neighbour also. 
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CHAPTSE III 


^/.M-c£ri07VX. £^xyp£jd ^£/m},cxirzciiJj:dofi3 

1 L - . Iri'croduclioa 

In this chapter we shall be dealing with structures 
oi terrahedral symmetry. Fost of the II-IV and III-V semi- 
conductors possess this symmetry. The electronic structure 
of semiconductors with impurities in substitutional and 
intei^titial position nas many unresolved aspects and is 
the focus of attention of many groups of scientists. The trans- 
lational symufitry of the crystal is broken by the impurity and 
the defect levels deep in the band gap of the semiconductors 
are difficult to investigate by ab-initio analytical theory, 
though they are important in understanding the electronic 
properties of the semiconductors. Hence, earlier attempts 
to apply this approach have been restricted to semi-empirical 
methods, specially the Extended Htickel Method, which is a 
nonself consistent procedure. It was proposed by Ifessmer and 
Vi/atkins [1] and Larkin [2 ]to surround the impurity by a 
cluster of substrate atoms and apply methods from molecular 
theory. One might expect that such a description could 
handle deep defect levels in semiconductors, which have 
localized orbitals and states associated with extended 

impurity complexes. Usually, levels in the band gap with 

E 

the position > from any one of the band edges, valence 

or conduction, is called a deep level, A level close to the 
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conduction band edge is called a donor level and a 

level close to the valence band edge is called an acceptor 
level. “8^ is the band gap. 

The cluster model is known to account well the electronic 
structure oi substitutional or interstitial impurities present 
in the solid. Since Messmer, Watkins and Larkin many authors 
have demonstrated the effectiveness of the cluster approach to 
deep impurities in semiconductors. Of the various methods 
applied, ilS Xq, SCF method stands out. 

In this chapter we will first tabulate the symmetry 
formed by the ligands and then go over to the FS Xa SOP 
calculations for systems like Si and InP, The impurities 
in Si have been extensively studied by other authors [3-7J. 

In our work the transition metal impurities in InP are 
considered along with the different charge states of Pe in 
InP, The effect of alkalirmetal impurities in InP is also 
studied, 

III, 2 Symmetry and Basis Functions 

th 

The elements of the group are in a special 

situation, since their electropositivity and electro- 
negativity are both 4. This comes from the configuration 

3 

of their outer shell electrons (ns) (np ), All elements 

except for Pb crystallize into diamond structures 

(Figure 3,1a), in which four electrons per atom tain part in 
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the electron pair bonds 5 and each atom extends four chemical 
bonds in the tetrahedral directions forming a three diixen- 
sional networks Such crystals are named as covalent crystals 
or nomopola,r-bonded crystal or valence bond crystal.. 

Ine one s and the three p orbitals can be linearly 
combined 5 if the energy difference between these orbitals is 
small, to obtain rS"! 

L 


T 
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1 

2 
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" Px 
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h 
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Pz^ 


These fuuctions have a large amplitude in the direction 
denoted by the subscripts. Prom these tetrahedral hybrid 
orbitals one can construct bonding orbitals and anti-bonding 
orbitals by linear combination with similar hyorid orbitals 
of the neighbouring atoms. When a pair of electrons with 
opposite spins is accommodated in a bonding orbital, a 
covalent bond is completed. 

The above feature may also be described in terms of 
band theory. Por a large lattice parameter the s and p 
orbitals of the valence electrons in all atoms are degenerate, 
at the level of an isolated atom. With decreasing lattice 
parameter, the s and p level each form an energy bands the s and 
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p bands can accommodate 2 and 6 electrons per atom, respecti- 
vely. lor a srriallcr lattice constant, the s-band 8.nd p-band 
cross each other and tney mix into two bands of completely 

different character. The lower band cor re- spends to the 

3 

baiid of bonding orbitals of Sjj hybridization and the higher 
to that of the anti-bonding orbitals (Figure 3.1b). Each of 
tn^ni contains four electrons per atom. 

Oonsidering only the first nearest neighbours, we now 
construct the symiretry i?ith respect to the center of a 
cluster or molecule. The irreducible representation of the 
tetrahedral symmetry group is indiceted on the left-most 
column of Table 3.1. The second column shov's the o-'ntral 
atom orbitals which form the basis of this representation, 

Tnird and fourth columns show the c- and Tc-combinations of 
tb) ligand atoms which have the same symmetry as the central 
atom orbital in the saire row. 

III. 3 Silicon 

The chairacte ristic of Si is its strong cova,lency with 
a large pile up of charge in the bonds between atoms. 

Tnis raises two essential questions on the applicability of 
iiS Xj., 3CF iiEthod to the molecular description within the 
cluster: whether the muffin-tin approxin at ion of the potential 

xs sufficient for the representation of real strongly non- 
spharical potential, and whether the local potential 



101 



.Fig.;3‘1b. ■ 

EnttftV bond of diomond stfuctufi as o 
f ynjctlon of jo td^tlce consttJOt . 



TABJE 3 . 


' I M^al 

Representation | Orbital 


Orbitals for Tetrah* 
Ligand o 
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approxiifiation of exchange and correlation is appropriate 
for a system with such large density fluctuations. Howeyer, 
this method has successfully been applied to a, number of 
molecules x\dth covalent bonds [9]o lo improve upon the charge 
density to a large extent, spheres without nuclear charge 
can be introduced in the interstitial position, where the 
potential is spherically averaged instead of assigning a 
constant potential to this region [10,11]. The local potential 
approxinia.t ion of exchange and correlation in Si as compa,red 
vdth screened excha,nge shows an error of approximately 15% in 
the self- energy while the error in exchange gap ennancement 
is 2-3% [12]. 

Our calculated one electron energy levels for Si, by 
Its SOP met nod ai’e indicated in Figure 3.2. The levels are 
filled according to the Pauli exclusion principle and are ..les-i-g” 
nated by the irreducible representation of the tetrahedral 
point .groupo The 2t2--level is the last filled level and 
3a.^ is the first empty state corresponding to tht; valence 
band and the conduction band edges, n' spectively. The 
energy difference between these two b^nd edges, called the band 
gap, is found to be approximately 2.0 eV which is in fair 
agreement with tkt’ experimental value of 1 .H eV [13]» 

We tried to look at the bend gap by changing the exchange 
parameter a. The value 0.72751 for is taken from Schwarts 
[1 4], The next two values were obtained by increasing a in 




relolive scol 



steps of 0.1 . The band gap is found to increase and the 
levels in the valence band goes deeper as the value of cc 
is increased. The charge distribution for different a 
in the different regions is indicated in Table 3.2. The 
corresponding one electron spectra is presented in Figure 3.3. 

In the process of calculation of Cartling et.al. [3 ],four 
interstitial spheres were introcuced in addition to the 
fi^/e Si spheres. There were 12 hydrogen atoms to saturate 
the dangling bonds. In the present calculation neutral Si 
atoms are put together to form the cluster and hence the 
charge on the Watson sphere was set at zero. Table 3.2 points that 


TABLE 3.2 

Cha rge Distribution in Si for the Different Exchange 

Parameters 


Eegion 

< a=0. 72751 | 
1 } 

a =0,82751 

1 a;=o. 92751 

1 

J Cartling [3] 

Extraraole cular 

2.5876 

1 ,2557 

1 .787* 

0,19 

Central ion 

12.9365 

12,4702 

12.9573 

13.2 

ligand ion 

11 .6007 

13.7056 

12.6590 

12.3 

Inter sphere 

8.0695 

1 .4505 

4.6163 

8.39 
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the charge on the central ion is nearly the same. This 
shows that the central ion charge is not affected much 
for a slight change in the boundary condition 5 after the 
self consistency has been achieved. 

A vacancy is then created by removing the central Si 

atom. It is seen that the t 2 “level moves into the band gap, 

of 

This is because of breaking^onds betv/een the central atomic 
site and its neighbours. So this 2 t 2 “level in the band 
gap of Si spectrum is associated with the broken bonds at 
the vacancy site. An impurity may be inserted in the vacancy 
to study the other defects in the 3 am.ple. A simple argument 
for why the a-j-state ought to be lower than the t 2 “State [15] 
is illustrated in Figure 5,4, Here the atomic valence s and 
p states are considered to interact with the a^ - and t 2 '"States 
of the lattice vacancy into x^rhich the impurity is inserted [5], 
In the present calculation a shallow impurity level in the 
gap at + 0.55 eV is obtained. The results are shown both 
in absolute and relative scale in Figure 3.2. In the 
relative scale the upper most completely filled levels are 
matched together. The 1 a^~ and 2a^- levels move up towards 
the valence band edge with the e-level forming the valence 
band edge. In this five atom cluster no effect of dangling 


bonds was considered. 
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Schematic representation of vacancy-impurityCcnd 
vacancy-silicon interactions where(a) represents 
crystalline silicphj[b)the lattice vacancy, and(c) the 
substitutional impurity. 



in Indium Phosphide 


The electronic structure of transition metal (TM) 
impurities in semiconductors are of much interest both from 
ph^^sical and technological points of view. These dopants 
are introduced deliberately as charge compensators. The 3d 
transition metal elements form localized states, called the 
deep impurity state, in the band gap of the host semiconductor. 
These states play an important role in modifyiiig the 
electronic properties of these semiconductors, For example, 

Or makes GaAs semi-insulating and Fe plays a similar role 
in InP. Or doped GaAs is used as semi-insulating substrate 
for epitaxial growths and also in devices like field effect 
transistors and photoconductors. The TM doped InP have a 
very promising technological future since the devices made 
out of them are much faster than those made out of Si. 

A battery of experimental tools such as electron 
par-amagnetic resonance (EPR) , luminescence, transient 
capacitance etc. are employed to study the TlVl impurities in 
semi-conductors [16] , But the results are difficult to 
analyze if no accurate energy level structures are available. 
The effective mass theory rightly accounts for the shallow 
impurity levels in the band gap but there is no such theory 
to take care of deep levels and the detailed understanding 
of TM impurity states has continued to elude us. 



In this part of the thesis we discuss the electronic 
stiuctures of the neutral, unrelaxed, transition metal 
impurities liJr® Or, i.n, Fe, Co, Ni and Cu in Inf, It is 
well known from EPR experinents [16] that the T1 atom 
occupies the cation site (In-site), due to chemical cons- 
traints, in III— V semicond 'lectors w'ith a minimal of lattice 
distort ion» The case of single In va,cs-ncy in In? is also 
t ake n up » 

The exchange parameters used in the 5/iS SCF calculations 
are due to Schwartz [14] « Table 3»3 gives the cc -values of 
the host atoms and the impurity atoms® 


TABES 3®3 

u -Values of the Atoms Used in Our Calculation [141 


Atom 

T 

! 

1 a 

t 

f 

In 

0.70000 

p 

0,72620 

Cr 

0.71352 

Ml 

0,71279 

Fe 

0.71 151 

Co 

0.71018 

Ri 

0.70897 

Gu 

0.70697 



The bond length of InP is 4.08 a.u. [17] and so the 
central ion, ligand ion and the Watson sphere radii were 
fixed at 3.164 (muffin-tin radius of In [18]), 1.916 and 
6.996 a.u., respectively. 

To start with^we consider the 'perfect cluster' where 

the In is placed at the center of the cluster with four 

neighbouring P atoms placed in a tetrahedral configuration. 

Then a vacancy is created by removing the central In atom 

to 

but v/ithout letting the P atoms^^relax . The TM eleuBnts 
are put one after the other to study the electronic 
structure of the 'defect cluster'. There are 3 valence 
electrons available from In and 5 from each of the P atoms. 
Since each of the four In-P bonds uses two electrons of 
opposite spins, a net charge of (3 +4x5 -8=) 15 electrons 
is pushed up to the Watson sphere to maintain the electrical 
neutrality of the cluster. The uppermost filled level and 
the lowerfli'ost empty level correspond to the valence and 
conduction band edges. Since the bonding of the P atoms 
of the cluster with the second nearest neighbour remain 
unchanged for any substitutional impurity considered within 
the cluster, the charge on the Watson sphere was kept 
constant to simulate the solid. 

The calculated energy level spectra of the clusters 
are shown in figures 3.5 and 3.6, respectively. The 

tetrahedral symmetry and hence the 


clusters possess a 




otf'^ energy in cV for the different transition metal 



OiH 



Relative energy in eV for different trar^sition metal impurities in InP. 



energy levels are designated by the irreducible represen- 
tation of the tetrahedral point group The energy levels 

are filled up in accordance with the Pauli exclusion 
principle. Pigure 3.5(a) correspond to InP where the last 
lilled level is of t 2 -syriiKietry and the first empty level 
is of a-| -symmetr ; . This symnntry is also found in the 
bulk band structure calculations where the top of the valence 
band is p-like and the bottom of the conduction band is 
s-like « The difference in energy between these two levels, 
called the band gap, is 1,08 e? and is in reasonably good 

agreement with the experimental value of 1.35 eY [13]. 

the 

It is important to note that/f/jS X^; SOP method, in which no 
empirical parameters are employed, is able to reproduce 
features of the bulk band structure. Similar accuracy has 
been reported by Hemstreet [ 1 ], Pazzio et.al. p9] 

De Leo et al. [5] . Tnis may be contrasted with the molecular 
orbital approaches [20-21], where the band gap is over 
estimated by atleast a factor of two inspite of empirical 
p arame t rizat ion . 

We now create a vacancy at the In-site. Pigure 3.5(b) 
shows that the tg-level moves up into the band gap from 
the valence band edge. This t 2 “state is usually designated 
by the vacancy dangling bond hybrid. This six-fold 

degenerate contains just three electrons as the 

three valence electrons lost by the removal of central In 
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atom. This t 2 "level is an acceptor level in the gap 
associated with the broken bonds due to creation of vacancy. 

The energy diiference between the t^^^ level and the valence 
bahid edge is 0.4 eV. Theoretical calculations of vacancy 
in OLher semiconductors predict a partially occupied level 

in the gap [5 519 ]. 

The vacancy created by the removal of In is now filled 

by Or. The energy levels are shov/n in Figure 3.5(c), The 

atomic 3d-orbital of tte TM is split due to the tetratedral 

environment into two crystal field resonance (CFR) levels, 

t 2 and e . This bonding t 2 has an anti-bonding 

counterpa.rt in the form of t^ . At the high Z-limit , 

\\fhere the atomic 3d levels lie far deeper than the anion 

derived host t 2 oi'bitals, there is little interaction between 

DBH 

them. Hence the p-like anti-bonding tg orbital is Qust 
above valence band maximum. We can see the appearance of 
these CFR levels in Or doped InP (Figure 3.5(c)). Four of 
the six electrons of the Or impurity are used up in 
healing the broken bonds created due to vacancy. The 
bonding orbital is now expected to acquire a s-d character 
in contrast to the s-p bonding in InP. But this is not 
very clear from our calculations. The partially filled 

moves into the band gap while tj is at the valence 
band edge. The unoccupied tg^^ level is pushed up almost 
into the conduction band. The e— level in the gap is at 
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+ 0.41 eV and the crystal field splitting is 0.68 eV 
( d ifieirerice between partially filled e-level and the 
unoccupied t 2 “level). The level is just 0.14 eV 

below the unoccupied 2a-)- level at the bottom of the 
conduction band. 

CT^R 

When Gr is replaced by iJn, the unfilled e " levels get 
filled and is pulled down to form the valence band edge 


(Figure 3.5(d)). The empty t™ level is at + 0.65 eV. 


"Y 


Ivmi is then replaced by Fe. Figure 3.6(a) shows the 


1 ai 


and tte t?^^ levels to be pulled further down in the 


,I)BH 


valence band compared to that of B/n. The t 2 level in 

CFH ^ +.DBH 

the band gap contains just one electron. The e and. ^2 
energy difference is 0.54 eV, and the t 2 2a-j- are 

separated by 0.54 eV, So the t™ level lies exactly in 
the middle of the band gap and this provides a clue to 
the semi-insulating nature of Fe in InP [22]. 

The replacement of Fe by Co, Ni and Cu shows a monotonic 
trend in the energy level structure. Comparing Figures 
3.6 (a,b,c,d) one finds that the energy difference between 


is 


almost same 

level increases 


CRR 

the two levels e - and 2a-)- 
in all the cases. The occupation of t 
from one to four as one goes from Fe to Cu. The 1a^- 
ard t™ levels go deeper as the atomic number Z increases 

due to the InoreaelBgly attractive potential of the heavier 
TM impurities. The other notabile feature is the gradual 
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in position of the level in the band gap. The 

energy difference between the and t^®^ level is 

into eV for Co while that for Fi and Cu are 0.41 and 0.30 eV, 
respectively. Tnis shows that Cu is a shallow acceptor 
\,naile fe is a deep level. In the relative energy scale 

fcqin 

tne j-'osition of the e level is found to be fixed « The 

CIPl 

e level forms a a~bond with the second nearest neigh-- 

bour« Our cluster does not have In as the second nearest 

neighbours. So tnere is no stabilizing mechanism in our 

Qrp 

calculation for e level,, Hence it does not get pulled 
down in energy as opposed to a-j- and tg levels. 

A brief survey of other cluster calculations is ''.'orth- 
while at this juncture. The study of TK impurities in 
Ci [7 j 2'3]> GaAs and GaP [25] shov/s a similar trend 

as above the level in the gap, 

Pe leo et.al. [5] have sho^vn that the electronic 
structures appropriate of the Si 16 cluster are closely 
similar to those of the smaller Si 4 cluster. In fact 
there is very little change in a-^ - and 12 “ state positions 
relative to each other and valance band maxima in going from 
Si 4 to Si 16. A first principle calculation of Pe in InP [22] 
also predicts the iampurity level (S,^ + 0.7 el) to be at the 

mid j-ap (0 = 1 «4 eV) which is consistent wath our calculation 

S 

(E + 0,54 eV with E =1.08 eV). This proves that cluster 
truncation effects do not alter tte electronic structure of 



I'e ir. in.;-' signii leant ly. It provides confidence in our 
c a i c ul at ions. 

■ igux'o j /J shows the cha.rge residing on the central 
iori, as obtained by L3 calculation for any energy 

' O "R^'R "R ~Pi 

level ia irnasure of localization) ^ for 0 “"^ and t 2 states, 
hion the localisation is plotted against the atomic number, 
it secorr.es evident that the impurity levels of Cr have more 
atomic character than the other TM impurities. For example, 
tnet.y" level in Cr has a localization of 40f as contrasted 

Ti'Drr 

to i,’ tor Cu. The reason is that for Cu, is valence 

oaiid lii.e and nance the wave function may be expanded in 

tonus of valence Bloch functions using effective mass 

t;r.Oi\,'. Tntorestingly, Fe marks the transition point. The 

DBH 

iocai i.-at uon of the wave function for t 2 shows an abrupt 
v.ariation at .'Je. 

■inally, Figure 3.8 shows the charge transferred from 
the central ion region to the ligands and/or into the 
intors'.'here region. The charge lost by the TL: goes mainly to 
the intersphere region (30 ''j) and less to the nearest 
neighbour F atoms (20%). The charge transferred to the P 
atoms from re, Co, Fi is constant (-17.5^). This is in 

corjj'orrr.ity with electronegativity (X) consideration. Fe , 

Co and Ki have the same electronegativity (x = 1 .8) while 
P has a higher electronegativity (X = 2.1). Similarly, 
ivm which has X= 1.5 transfers a larger amount (25%) of its 
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Fig. 3 7 

voriotlon of loeolliotlon for th* differcM trotrsitlon 
motol impurUlts lt> InP In orWrory unitt. 
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Fig. 3*8 

Charge Trcnsfdf from the transition metal impurities 
in tnP. 
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chai'ge to P„ This shows the itinerant character of the 
non-core electrons and it is interesting that even a 
siniple cluster as ours oould hring this out so well. 

Similar trends are shown by other Th. impurities [22 

The charge transfer does not show a monotonic behaviour. 
It is small for Fe but larger for two neignbouring elements 
tin and Co, This is certainly in conformity with the 
semi- insulating behaviour of Fe in InP, Our energy level 
calculations corroborate with this fact. The tg ^ level 
is in the mid-gap for Fe thus exhibiting extreme localization. 
The same t^®^ level is closer to conduction band edge for 
j.n and to the valence band edge for Co, They thus hybridize 
v’ith the host and lose sons of their charge. 

In recent years careful electrical studies on Fe and 

doped IrJ? have been conducted. The resistivity of Fe doped 

g" 3 5 

InP is > ^0 ^ range 10-10 0 cm 

for Go doped In.?, This we attribute to the fact that the 

level is in the mid-gap and hence highly localized 

for Fe while it is closer to the valence band edge for Co, 

Indeed, it is experinentally found that the main deep 

level is at + 0.75 eV for Fe [26] , while it is at 

^ + 0,32 eV [ 27 ] for Co. Experiments show one more level 

in the band gap, Ttese are not interpreted by our theory 

and we feel that these may be associated with the different 

charge states of Pe 3Xid Go® 
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-fe conclude by stating that we have successfully 
identified the main deep levels of Pe and Co doped InP 
and interpreted the basic features of their electronic 
structure and properties. 


III«5 -different Charge States of Pe in InP 

After finding out the electronic structure of TM 
impurities and comparing our results for Pe with other 
theoretical [22] and experimental results [26] , we conclude 
that Pe makes InP semi- insulating. The obvious question 
arises v/nat happens when Pe is present as an impurity in 
Ini in different charge states? hfe shall now study the 
el'vctronic level structure for the different charge states 
of Pe . 


3 

Th--^ central ion is linked with the ligands by sp 
h^z-bridizat ion in these semiconductors. The different 
chargi- states correspond to different number of electrons 


in the outermost 


shell of the atom or ion. Keeping 


the sp"" hybridization in mind, we designate the different 
charge states of Pe in InP as Pe^'^ (3d^4s^), P9^'^(3d^4s^ ) , 
and Pe 




c 2 

(3d 4s ), respectively, e 




Pe(3d^4s^) = Pe(3d\s4p^) == Pe(3d^ = Pe^" 

Pe(3d^4s^) = Pe(3d^4s4p^) = Pe(3d^ = Pe^ 
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1 Xl-h 

j. he liOtatiDii yd merely implies that there are 6-n 3d 
wave lUiiCcioriS centered on j?e , Similar nomenclature has 
been used lor different charge states of Or in C-aAs [28^29], 

i he theoretical one electron spectrum of the 

diiiurent charge states of le in InP, calculated by us by 

IS IciddH' method, is shown in Figure 3.9o The 3d level of 

C'R'R 

?e breens up into e "■ and t^” levels. The level 

forms the valence band edge for all the three charge states of 

_nT5H 

7e in xnP. The t^' level, arising due to the broken bonds, 
in the hand gap contains 2,1 and 0 electrons for Fe 
Fe , and J'g , respectively. This level moves towards the 
^'aionce fand edf;e with increase in population and its 
;>os:Ltion in tie gap is at i 0„97j + 0,54 and 

■' O.F'j eV, respectively (Figure 3,9), 

The normalized population density analysis (Figure 

3.1'’)) shovv's the localization, in arbitrary'’ units, of tg^^ 

C '<’7 

and e " levels. This trend is similar to that shown by 
other Tii elements used as dopants in III-V semiconductors 

[7,19,24,25] . 

The semi "insulating behaviour for neutral Fe in 
InP has already been pointed out [22,30]. The study of 
the charge distribution in the different regions of the 
cluster (Table 3.4), shows a decrease in the central ion 
charge in going from Fe^^ to Fe^^ or Fe^"^. The ligand charge 
is found to remain almost constant. The charge transfer 
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Vorlotion of ioccrizotion for the different 
chorge states -of Fe in InP if^ orbitrofy 
units. 
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ix'orr. tho central le atom to the nearest P atom is approxi- 
matel/ , constant in all the three cases. The 

extra charge goes to the intersphere region. 


TABIT 3,4 

Charge listrlbution for the Different Charge 
States of Pe in InP 


Region 

; Pe^ ^ 

t 

f 

9 - 

1 Pe ' ' 

I 1 

j i 

Pe^+ 

Ixt r amo le cular 

0,044 

0,017 

0.036 

Central ion 

20,908 

24,234 

20,063 

ligand ion 

11 .050 

11 ,103 

1 1 .530 

liitersphere 

4,648 

2,334 

3.778 


This proves that 
or re 


^ 2 + . 

localization in re is 


more than Pe 


+ 


In addition to these, the charge density at the 
origin was also determined and used to calculate the isomer 
shifts of the charge states Pe'^ ‘ and Pe^ of Pe in 

InP. As yet no Fdsshauer data of Pe in InP is available. 
Table 3.5 shows the isomer shift of different charge 
states of Pe in InP calculated by us along with the experi- 


:n vvf* o 1 


n-P -1 c? r^TTTO -r* at 


'XTO'ni ]nAc?i“o 
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TABLE 3.5 

Isomer Shift of Be in Various Hosts 


Theoretical Experimental 

J 

(present calculation) ' 

- - t f 


( in 

In?) 

0.63 

Go 

0.04 

Ge 

0.331 

(in 

InP) 

0.37 

Ki 

0.02 

In 

0,410 

(in 

InP) 

0.62 

Cu 

0.225 

Si 

0.505 


The tteoretical as v/ell as the experimental values are 
relative to metallic iron. The general trend is that for 
metp-llic hosts which are in the neighbourhood of Be in 
the ])eriodic table, the isomer shift is small compared to 
tho case of semiconductors like Ge , In and Si. This is 
because the charge density distribution for Co and Ni should 
almost be identical to that of Fe because of their proximity 
in the periodic table. The semi-insulating nature of InP 
has cxli'eady been pointed out with the —level lying in 
the middle of the band gap. The charge on the central ion 
for (Table 3.4) shows that Fe retains more of its 

original atomic character, fe know from neutron diffraction 
experiments tha.t in the metallic iron, the iron charge 
distribution is substantially similar to its atomic counter- 
part [31 ]. Hence the isomer shift of Fe in InP with 
respect to metallic iron would be expected to be small as 
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compared to Fe in In and Si„ Por the other two charge states, 

1 + 3 + D'RH 

Pe and Pe , the t' 2 is either closer to the conduction 
band or the valence band edge. The wave functions in these 
two cases can then be expanded in terms of conduction band 
and valence band Bloch states. These Bloch states have 
mostly In? cnaracter and only maiginally metallic iron 
character. Hence the isomer shift is large. We can look at 
it somewhat differently. A la.rge amount of charge goes out 
of the central-ion (see Table 5o4) thus decreasing their 
contribution of charge density at tte origin. Hence the 
isomer shift increases depending on the charge given 
out . 

The low energy photoconductivity experiments [32] 

1 + 2 + 

find the Pe to Pe transitions at 0.34 e¥ which is in 

excellent agreement with our calculated result of 0,36 eV. 

2+ 3 + 

However, the Pe to Pe transition has an experimental 
value of 0,64 eY which is in contrast to our value of 
0,33 eV, We expect that the detailed discrepancies may be 
removed by expanding the cluster to the second nearest 
neighbour of the impurity. An increase in the size of the 
cluster will no doubt change the band gap [5,33], but the 
general trend of the electronic levels is expected to remain 
identical. However, experiments like Hall effect and 
deep level transient spectroscopy [27] indicate additional 



deep levels both in Fe and Co doped InP when the samples 
are subjected to heat treatment at 350‘^C for 1 minute, Vfe 
leel that only a spin unrestricted calculation can lead to 
more than one state in the gap. 

To the best of our knov/ledgej this is the first 
attempt made to understand the detailed electronic structure 
of the different charge states of Pe in In?, The above 
general agreement with the experimental data seem to be 
promising. We expect that calculated isomer shift results 
will stimulate further experimental activity in this area, 

III. 6 Electronic Structures of Substitutional Alkali-fetals 

in InP 

The interaction of the alkali-metal impurities with 
defect dangling bonds in semiconductors is a sub;]ect of 
current interest and controversy. The importance of under- 
standing the hydrogen vacancy interaction is clearly suggested 
by simple chemical reasoning: The vacancy or other defect 
dangling bonds, associated v’ith defect electrical activity, 
could be passivated by appropriately located hydrogen 
impurities. In fact there is experimental evidence for the 
neutralization by atomic hydrogen dangling bond related 
electrical activity in Si [34]. This effect has importance in 
solar cell and electronic device technology. Such chemical 
reasoning and observations therefore suggest the following 
question: Can other group-I impurities with their characteristic 



single unpaired s-ele ctrons ,, also passivate defect-related 
danpling bonds in other semiconductors'’ 

Alkali-iriet als are normally observed in semiconductors 
at the tetrahedral interstitial site [35];. and when located 
there it is well understood. In irradiated material, however, 
impurities could be trapped at lattice vacancies. Therefore, 
one might expect to find a variety of alkali-metal- vacancy 
complexes, Experiments suggest the existence of alkali-metal- 
vacancy and other alkali-metal defect in complexes of Si[34]. 

As yet no experiTCntal data for alkali-nBtal- vacancy complexes 
in InP exists, 

Vfe proceed here to explore these questions of electrica 
activity and dangling bond passivation by calculating the 
single particle electronic structures of substitutional li, Fa, 
and K in In?, Tetrahedral symmetry is preserved in all cases 
in these calculations. 

According to the model proposed by Bernholo et,al,[36], 
dealing with Group-I substitutional i.mpurities in Si, the 
group-I valence s-states tend to repair the broken a^ - 
symmetry dangling-bond combinations associated with the 
lattice vacancy. Lacking a 'p' valence electron, however, 
the t 2 electronic structure of the vacancy will tend to 
remain int a ct , 

The electronic structures for the substitutional 
alkali impurities, Li, Fa, and K, shown in Figure 3,11, have 
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Relative CTttrgy in ev tor pure lnP,Q vacancy at me in-siie an 

alkoU -metal impufities Ifke Li>No ond K tn InP. The filled' circles 
denote the number of electrons in the particular state. 



132 


certain similarities, i'e find here a vacancy like tg-state 
in the gap which is occupied by four electrons. This, 
tnerefore^is similar to the case of a vacancy (Pig,3.1Hb)) 
at the In-site. The extra electron in the tg^ga-P state, as 
compared to the vacancy case, is provided by the s electron 
01 the aikali-me tals . The 1a^- state forms the valence band 
edge as in the vacancy case. This is in contrast to 
the la-j -state in pure InP which is Bloch like. The position 
of t2-gap levels are at + 0.4, + 0.29, + 0.14, and 

+ 0.71 eV, for a vacancy at the In site, and respectively 
\\/ith Li, Ba, K. as impurities. 

Be leo et al. [5] proposed a simple model for the 
alkaii-inetal impurity in Si , showing the changf/shifts in the 
a.j - and t2- levels (Figure 3.12). The computational results 
suggest tha.t there is only a relatively weak interaction 
between the atomic valence ns and np states. These inter- 
actions are, therefore, similar tothoee of the vacancy- 
dangling-bond. The alkali-qietal core states being substantial! 
deeper, hardly play any role in the alkali-metal- vacancy 
interaction. Physically, we say that the defect state has 
been 'pinned' in the gap. If the impurity had introduced a 
near band gap 'p' state, then we would have expected to find 
the gap swept clean of the to states as the impurity 'p’ and 
vacancy tg-states interacted [37,38]. Instead, we find that 
the vacancy t2 state remains almost intact® 
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VACANCY SUBSTiTUT!ONAL ATOMIC 

ALKALI ALKAL! 

Fig. 3*12 


LiQcind-ficld model for th®8l®ctror4ic structure of the 
substitutional alkali-metal impurUiesintermsof the 
interactions between the vacancy orbitals and the 
alkalf-metal-atom orbitals. 
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Table 3«6 shows the charge distribution in the 
different regionsof the cluster^ The charge on the ligands 
IS almost constant for the different alkali-metals as 
impurities. The charges on li and Na indicate that they tend 
to retain their atomic character. A major part of the 
cnarge transferred from K goes to the intersphere region. 

The intersphere cnarge increases as one goes from li to K, 
through Ka, 


TABLE 3»6 


The Total 'Electronic Charge in the Different Regions 
o f the Cluster with Alkali-ketals as Impurities 


I ! I 


Region \ 

? 

f 

Li ; Inf 

! Na : InP 

! I 

J - - -? 

K;InP 

Extra molecular 

0.036 

0,037 

0.028 

Central ion 

3.038 

10,173 

15.099 

Ligand ion 

11,165 

11 .160 

11 .461 

Inter sphere 

0,226 

1 .0Q8 

3.029 


The electronic structure and the charge distribution 
suggest that these impurities x;ould not passivate :he vacancy 
dangling bonds. In fact, the impurity system remnin electri- 
cally active as a multiple donor. 



135 


REFERENCES 


1. R.P.Nessmer and G.I) jA/atkins , 1970, Phys » Rev, Lett , 656 , 

2. E«P, Larkin, 1971, L. Phys. C 4, 3065, 3077« 

3. B»Gartling, B.Ross and U.Wahlgren, 1 973, Chem, Phys , Lett, , 
il, 380, 

4. LoCartling, 1975 , Phys, 08, 31?, 5183, 

5. Gary G.De Leo, N.B. Bowler and G.L. Watkins, 1984, Phys. 

Rev. B 29, 3193, 181 9, 

6. C-ary G, Be Leo, G.D, Watkins, V/.B. Bowler, 1 982 Phys, Rev. 

B 25, 4962, 4972. 

7. L. A.Hemstreet , 1977, Phys, Rev, B. 1^, 834, 

b, T. Mat Sahara, 1982, The Structuxe and Properties of 
Matter ( Springer- Verlag Berlin, Heidelberg, 

9, J . '>'.1 , Connoly, H, Siegbahn, U, Bellas and C.Nordling, 1973, 

J, Chern. Phys, 4265, 

10. J, Keller, 1971, J, Phys. C 4, 185. 

11. A. Pl. W illiams and J.Van W.Iv organ, 1972, J. Phys. C L293. 

12. E.O.Kane, 1972, Phys. Rev. B5, 1493, 

13. C.Kittel, 1977, Introduction to Solid State Physics, 

5th Edition, p.2l0, Wiley Eastern Limited, New Lelhi, 

14. K.Schwartz, 1972, Phys, Rev, B 2466. 

15. G.B.Bachelet , G.B.Baraff and F.Schluter, 1981, Phys. Rev. 

B 24, 4736. 

16. U.Kaufmann and J, Schneider, 1980, Bestkorper Probleme , 

20, 8787. 

17. H.B.Wolf, 1971 , Semiconductors, Wiley Interscience, New York 

18. V.L.Moruzzi, J.B.Janak and A.R. Williams , 1978, Calculated 
Electronic Properties of iVAtals, New York, Pergamon, 

19 . A.Bazzio, J.R.liete, M.L.De Siqueira, 1979, J. Phys. 01^, 
513 , 3469» 



20« 


1 36 


C. Weigel, l.Peak, J.W. Corbett, GoD. Watkins, and k,P. 
Jfessmer, 1973, Phys. Rev, B 8, 2096, 

21, Vijay A. Singh, C. Weigel, J.W. Corbett and 1 . 17 , Roth, 1 977, 

Phy, Stat, Sol. (b) 637, 

22, 7^i;iay A. Singh and A.Zunger, 1933, Bull, Am. Phys. Soc., 

28 , 288, 

23, R.KLeinhenz, V'ijay A. Singh and J.W. Corbett 1979 
(unpublis hed ) . 

24, A.Bazzio and J.R.liete, 1980, Phys, Rev, B 4710, 

23, N, Gemma, 1984, J, Phys. C 1_7, 2333. 

26, S.Bung, R.J. Nicholas and R.A, Stradling , 1979, J* Phys, 

0 12, 5145. 

2'7. ivi. S . Skolnick, R.G. Humphreys, P.R. Tapster, B. Cockayne and 
W.R.MacEwan, 1983, J, Pnys. C 16, 7003. 

28. J.J. Krebs and G.H.Stauss, 1977, Phys. Rev» B 1_^, 17. 

2 j. 1, A.Hemstreet and JoOalimmock, 1979, PhysiRev. B 20, 

1527. 

30. P.K.Khowash, D.C.Khan and Vijay A. Singh (submitted). 

31. W. Marshall and S.W.lovesey, 1971, Theory of Therical 
Neutron Scattering, Oxford, Clarendon Press. 

32. L.Evas, A.W.Smdth, P.J. Williams, B. Cockayne and 
'.W R. MacEwam., 1981, J. Phys. C l£, 5036. 

33. J.R.liete, V.M.S. Gomes, l.V.C.Assali and l.K.R. Scolfaro, 
1984, ICDS (Private communication). 

34. J.l. Benton, C.J. Doherty, S.D. Perris, D.l.Plamm, 
l.C . Kimerling and H. J.Ieamy, 1980, Phys. lett. 

670 . 

35. J.W. Corbett and G.D. Watkins, 1971, Ra,diat ion Effects in 
Semiconductors, Gordon and Breach, tlY. 

36 « J.Bernholc, N.O.lipari, S.T.Pantelides and M, Scheif ler , 
1982, Phys, Rev. B. 26, 5706, 

37, HHP. Hialmarson, E.B. Allen, H.Buttner, and J.D.Dow, 1980, 
j. Yac. Sci, Technol. H, 993- 


TT-.' 


D ,-Y*1 


T\ J Wn1 f*n7^H J "Hovr . 1 QRO » 



OHAPTES IV 
^I^on£jtu6lofn, 

IV.1 Aim 

to 

The aim of the present work has teen^ study theoreti- 
cally the electronic structure of a number of ionic y 
covalent and semiconductor solids, using MS S^^SOP method. 

The valence electron energy levels and the equilibrium charge 
distribution of all these systems have been analyzed to get 
informations on their bonding characteristics, optical proper- 
ties and electrical conductivity. 

The MS XqSCP method considers a cluster scooped out 
of the solid and replaces the effect of the rest of the solid 
by putting proper charges on the boundary of the cluster. 

This is expected to simulate the localized physics of the 
solid to a good approximation. Hence the method is very 
suitable for the study of the bonding in the mixed ionic and 
covalent systems and the localization of the impurity levels : 
in semiconductoTB. 

IV.2 Achievements 

Ife applied to MSZ^SCP method of calculation to 
systene like.ionic solids (Nadi, liP., NaP), penovskites 
(EMP^ (M^s Pe,:. Ni)) , ..ibtilea- Go, .Ni, ..Za)), 

metallic oxides . AlChg 
wit h t he SOP energies t he . chargee within, .each at oi^d, n^gion 



are also calculated. A charge transfer diagram was defined 

to get an understanding of the cheniica.! bonding (ionic and/or 
covalent) in these systems. The 3d degenerate level in the 
cation breaks up into t 2 ^ and levels in the ligand field, 
ihis crystal field splitting, 10 Dq, v/as determined for 
all the crystals with 3d cation and compared v/ith experiment 
'A7he never data were available. The transitions between two of 
the valence levels were calculated whenever we could locate 
the corresponding XPS or optical spectroscopic experimental 
data, finally the ^ibssbauer isomer shifts were calculated 
for iron containing samples. 

We started with some ionic solids, e.g. fTaCl, Naf 
find Lit. The self consistent charge distribution study 
showed a transfer of charge from the central ion to the 
ligands along with some transfer to the interatomic region. 
This showed that these well known ionic compounds are not 
perfectly ionic but show a marginal covalent character, 
for haCl the calculated charge distributions are in very 
good agreement with the experimental results. The energy 
levels of lif and Naf have the characteristic that they 
bunch into two groups with predominantly I? 2s and f 2p nature. 
The energy difference of the P 2s~like and P 2p-like states 
in liP is 1 .43 Ry as compared to the experimental result 
of 1.54 Ry obtained by X~ray photoelectron spectroscopy. 
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In the next stage the energy levels and the charge 
distributions were calculated for the fluorides and the 
oxides of different transition metals* The agreement of 
the energy difference of the F 2s like and P 2p like states 
in Pep2»^»36 Ry ,vrith the XPS result 1*54 Ry was encouraging. 
The crystal field splitting, 10 Dq, has also been calculated 
for the transition metal fluorides. The Ilg 2s and the Mg 2p 
levels were identified and their energy difference 2.85 Ry 
was found to be in fair agreement with the XPS data 2.58 Ry. 

The Ifdssbauer isomer shifts in all the iron contain- 
ing systems namely KPeP^, Pep2 and Pe alloy were 

calculated and compared with the existing results from the 
experiments. 

The first step towards the analysis of the covalent 
nature of the systems was to calculate the covalency 
para,meter X from the coefiicient of the ligand term in 
the extended eigen-function of upper most valence level. 

These values were matched with the obta-ined from the 
analysis of the neutron diffraction and nuclear magnetic 
resonance experiiDental data in terms of the molecular orbital 
theory. We get satisfactory agreements for KPeP^ and KRiP^, 
Por KPeP3 our results are = 5.58fo as against the neutron 

diffraction result = 6.24;^. Por ^3 they axe 7.55% 
and 5.11?3 respectively. The value of for MgO is 6,6? 
percent which is in good agreement with the nPROR 
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va.lue of 6.50 percent. 

¥e then brought in the concept of pseudoatoms 
and calculated the charge transfer from and/or to these 
a,toms and plot the charge transfer diagrams for all the 
solids studies by us. With the help of these diagrams we 
defined the measures of ionic bonding strength ard covalent 
bonding strength for these systems, showed that the 

na.tural state of the non-metallic solid is a mixed ionic- 
covalent state. In KaCl the ionic is the dominant bonding 
and covalent bonding is marginal whereas in PeP2 covalent 

bond strength is very large and ionic bonding is marginal. 

I'he transition metal perovskites and rutiles generally 
fall within these two limits. However j they are predominantly 
covalent. The oxides of transition metals also show very 
large values of covalent, bond strength with different degrees 
of minor ionic character. I-IgO is an exception to this 
general character in that it has a small charge transfer 
from the anion to the cation. This type of charge transfer 
is typical of semiconductors and we have referred to it as 
inverse ionic character, MgO is known to behave as a 
semiconductor. Inverse ionic charactor along with large 
covalency was shown by all the semiconductors worked out 
by us, e,g. Si, InP and transition metal impurities in 


InP, 



In the second part of the thesis MS Xi^SCI method 
was applied to semiconductors, especially transition metal 
doped InP, with a view to study the localization of the 
levels within the band gap. The cluster model had been 
applied to the well known semiconductors, such as Si and 
■G-aAs . It was found that Or makes G-aAs semi- insulating. 

expected that Pe plays a similar role in InP and worked 
out the self consistent field energy levels and their eigen- 
functions. Our calculations did predict a highly localized 
impurity level for Pe in InP. 


Por the other transition metal impurities in InP 
wG showed that the impurity level in the band gap moves 
closer to either valence baud or conduction band edges. 

Those states thus hybridize with the host and become deloca- 
lized as compared to Pe in InP. Very recent experiments on 
Pe and Co in InP confirms a higher localization of the 
former leading to a higher resistivity. 


kfe also performed calculations on the different charge 
states of Pe in InP. We found that the tg®^ -levels arising 
due to broken bonds, is exactly in the middle of the band 


gap for neutral Pe , whereas it moves 
for the other two charge states. The 
minimum for neutral Pe . The isomer s 


closer to the band edges 
charge transfer is 
hift for different 


charge states were calculated but unfortunately no experimen- 
tal data are available to compare with. However, these result 
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also confirm that neutral Pe makes InP semiinsulating. 

The electronic structure of the alkali metals as 
impurities in InP were also calculated^ It was found out 
that they interact weakly with the host thus pinning the 
impurity level in the gap. 

It is piedicted that G-aAs may soon supercede Si in 
technological importance. Since the devices made out of 
InP have proved to be much faster than those of Si, it is 
expected that InP is going to be a keen competitor to 
GaAs as a device material. Our work forms the theoretical 
basis for the analysis of the conductivity properties of 
transition metal doped in InP. 

IVe 3 limitation of the Method 

The cluster modal is known to emphasize the 
localized physics. The long range effects may require 
clusters too large to be readily handled. Accurate values 
of bulk cohesive energies, work functions and many effects 
associated wit n the fermi surface are outside the scope of 
even a cluster as large as 50 atoms. Ideally a cluster study 
of a property should be supplemented with the investigation 
of the convergence of results as a function of cluster 
size. However, this makes the ana.lysis very unwieldy and 
may not be within the computer economics of a researcher. 



143 


IV. 4 Scope for Further Work 

The present work concentrates on bonding character- 
istics and covalency in metallic compounds and deep levels in 
semiconductors. Since taese are all localized physicsj 
the consideration of nea-rest neighbours is justified. 

However, in semiconductors, the size of the cluster may 
be inciBased, By this the essential physics is not expected 
to change but the numerical values of the band gap will be 
in closer agreement with the experimental results. Hydrogen- 
can be used, in the second nearest neighbour position, to 
saturate the dangling bonds, in semiconductors. 

The xBuf fin-tin form of potential can te improved by 
taking non- muffin-tin corrections as in DV 1 q.SCF' method. 

Apart from these, a spin polarized calculation, if 
performed, c?.n lead to the actual magnetic state of the 
system. The magnetic form factor f(k) can then be calculated 
by taking the Fourier transform of the difference of the two 
charge densities 

ik.r 

i(k) = / dv [P^(r) - (£) p (4.1) 

The present activities of our group includes a spin polarised 
calculation for FeF2 » Also efforts are being made to extend 
the cluster to its second neighbours. 
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IV. 5 Conclusions 

We hope that the present work has extended our knowledge 
of the electronic structure of covalent solids and semiconductors 
with impurities. Quantitative measurements for the bonding 
characteristics have been developed for covalent solids and 
energy levels of the valence electrons have been determined 
forming the theoretical basis for optical studies, for 
semiconductors, the levels in the band gap with localized states 
have been studies in detailj interpreting photoelectric and 
electric conductivity. 






